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ABSTRACT 

Busek,  Robert  Henry.   M.S.C.E.,  Purdue  University, 
June  1971.   PRACTICAL  OPTIMIZATION  OF  STEEL  HIGHWAY  BRIDGE 
BEAMS.   Major  Professors:   J.  T.  Gaunt  and  A.  D.  M.  Lewis. 

A  program  for  the  optimization  of  a  rolled  section 
highway  bridge  girder  is  presented.   The  girder  is  designed 
by  minimizing  either  the  weight  or  the  cost.   The  cost 
function  used  is  a  representative  model  of  the  total  cost 
of  the  girder,  including  both  material  and  fabrication  costs 
The  function  presented  is  a  sample  and  can  be  altered  if  a 
different  cost  function  is  more  appropriate. 

The  program  produces  the  optimum  design  of  a  rolled 
section  girder,  using  either  composite  or  noncomposite 
design  theory.  •  The  required  cover  plates  are  determined, 
using  both  static  and  fatigue  stress  conditions.   The 
fatigue  stresses  are  based  on  500,000  cycles  of  stress. 
If  a  design  is  composite,  7/8- in.  shear  connector  studs 
are  provided  in  the  positive  moment  region.   No  connectors 
are  placed  in  the  negative  moment  region.   Ml  design 
theory,  analysis  theory,  and  optimization  constraints  are 
based  on  the  1969  AASHO  bridge  specifications.   The  method 
of  influence  lines  is  used  in  the  analysis. 

The  exhaustive  search  technique  is  used  as  the  main 
method  of  optimization.   The  cover  plate  thicknesses  arc 


X  1 


determined  by  using  the  method  of  interval  halving,  and  the 
locations  of  the  cover  plate  splices  arc  determined  by  using 
dynamic  programming. 

Computer  solutions  for  two  published  design  examples 
were  obtained  and  compared  to  the  published  results. 


CHAPTER  I 
INTRODUCTION 

1  .  1    Structural  Optimization 
Optimization  is  the  process  of  obtaining  "the  best" 
solution  to  a  problem  based  on  a  given  criterion.   Struc- 
tural optimization,  therefore,  is  the  process  of  obtaining 
the  "optimum"  solution  to  a  structural  problem.   This 
optimum  is  obtained  by  using  analysis  and  design  theory 
and  an  optimum- seeking  method  suited  to  the  problem. 
Structural  optimization  has  only  recently  been 
approached  systematically.   One  of  the  reasons  for  this  is 
that  the  computer  is  a  relatively  new  tool  for  the  engineer. 
When  doing  an  optimization,  there  are  a  large  number  of 
calculations  which  must  be  completed.   Performance  of  these 
calculations  by  hand  would  be  very  time  consuming.   The 
computer  gives  the  engineer  a  quick  and  fairly  inexpensive 
method  of  calculating  quantities  and  making  decisions  based 
on  these  calculations. 

Another  reason  for  the  recent  prominence  of  structural 
optimization  is  the  rising  cost  of  construction.   Since  the 
costs  of  material  and  labor  have  risen  steadily  over  the 
past  years,  and  since  the  trend  for  the  future  seems  to  be 
in  the  same  direction,  the  engineer  mast  place  mere  emphasis 


on  the  economics  of  design.   By  using  optimization  proce- 
dures, the  designer  may  actually  reduce  construction  costs 
by  minimizing  the  total  cost  of  the  structure. 

The  formulation  of  the  optimization  problem  depends  on 
the  type  of  structure  being  designed.   For  any  problem,  the 
design  variables,  the  design  constraints  and  the  criterion 
for  the  objective  function)  of  the  optimization  must  be 
defined.   It  is  the  values  of  the  design  parameters  which 
are  desired  so  that  the  value  of  the  objective  function  is 
a  minimum  and  the  constraints  arc  not  violated. 

The  design  parameters  required  depend  on  the  type  of 
optimization  problem  and  on  the  theoretical  equations  used. 
The  specification  used  to  control  the  design  specifies  the 
necessary  constraints  for  the  problem.   In  most  problems, 
both  equality  and  inequality  constraints  are  used. 

The  objective  function  can  be  any  desirable  measure  of 
the  effectiveness  of  the  design.   Some  of  the  measures  which 
have  been  used  are  weight,  cost,  cost-  benefit  ratios  and 
reliabil  i  ty  .   for  the  highway  bridge  girder,  minimum  weight 
has  been  used  for  many  years,  because  of  the  relationship 
which  exists  between  weight  and  cost.   In  this  optimization 
problem,  however,  the  cost  function  used  is  more  realistic 
because  it  includes  the  fabrication  costs. 

The  general  optimization  problem  can  he  expressed  in 
the  following  mathematical  form: 


Find  the  variables  x.     i  =  1 ,  2 ,  • • • ,  m 
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The  variable  vector,  x,  must  be  determined  so  that  the 
equality  constraints  (1-quation  1.2)  and  the  inequality 
constraints  (Equations  1.5  and  1.4)  are  satisfied  and  the 
objective  function,  G  (x) ,  is  optimized.   The  G  function  can 
be  either  maximized  or  minimized  depending  on  the  type  of 
problem . 

The  optimum- seeking  methods  which  are  available  to  the 
structural  designer  are  numerous.   Many  of  the  methods  are 
general  in  nature  and  can  be  used  for  any  type  of  optimiza- 
tion.  There  are,  however,  some  specific  methods  which  are 
best  suited  for  particular  types  or  formulations  of  problems 

1 . 2    Background  Review 
The  methods  of  optimization  can  be  divided  into  direct 
and  indirect  methods.   The  indirect  methods  are  mathematical 
methods  which  can  use  the  gradient  of  the  algebraic  form  of 
the  objective  function  to  solve  a  given  problem.   Direct 
methods  are  ones  which  make  trial  solutions  in  so:  e 
organized  manner.   1 ach  solution  is  compared  to  the  present 
optimum  and  a  decision  is  then  made  as  to  the  next  design 


step.   Each  of  the  types  has  its  own  advantages  and  dis- 
advantages when  used  in  a  particular  structural  optimization 
problem . 

Structural  optimization  problems  have  been  approached 
in  the  past  as  minimum  weight  design  problems.   Drucker  and 
Shield  [1  and  2]  have  presented  some  of  the  general  minimum 
weight  design  theory.   Faulkes  [3]  and  Ilcyman  [4]  have  used 
this  criterion  i-n  optimizing  frames,  and  Krishman  and  Shetty 
[5  and  6]  and  Haug  [7]  have  used  it  in  elastic  beam  design 
problems.   Ilahn  [8]  extended  the  work  to  include  plastic 
design  of  beams  and  simple  structures.   The  assumption  made 
in  these  minimum  weight  problems  is  that  the  cost  of  the 
girder  is  proportional  to  its  weight. 

Computer  design  methods  have  been  developed  for  highway 
bridge  design.   Sturman,  Albertson,  Cornell  and  Roesset  [9] 
developed  a  program  called  BRIDGE  which  can  be  used  to  make 
comparative  studies  of  different  types  of  bridges.   The 
system  can  also  completely  design  a  particular  structure,  hut 
it  makes  no  attempt  to  optimize  the  design.   Since  the  com- 
plete optimization  of  a  highway  bridge  is  a  long  and  compli- 
cated problem,  the  approach  so  far  has  been  to  optimize  the 
individual  parts  of  the  bridge  system. 

The  optimization  of  girders  lias  also  been  approached  by 
using  the  minimum  weight  method.   bdward  [10]  and  Holt  and 
Heithecker  [11]  have  attempted  to  optimize  the  girder  by 
minimizing  the  cross -sectional  area.   If  i  s  in  turn  minimizes 
the  weight  of  the  structure. 


The  total  cost  optimization  has  been  developed  for  both 
building  and  highway  girders.   Annamalai  [12]  uses  the  back- 
track programming  method  for  welded  plate  girders,  developed 
by  Lewis  [13],  to  optimize  the  design  of  building  plate 
girders  based  on  the  AISC  1969  specifications.   Okuba  [14] 
solves  the  girder  problem  by  using  a  linear  programming 
method . 

The  work  by  Razani  ,  Coble  and  DeSantis  [15,  16,  and  17] 
on  the  optimum  design  of  plate  girders  has  led  to  the 
development  of  a  computer-aided  design  system.   Coble  and 
DeSantis  [IS  and  19],  with  the  aid  of  the  Ohio  Department 
of  Transportation,  developed  a  program  referred  to  as  Girder 
Automated  Design  -  1  or  CAD- I.   The  program  completely 
designs  the  constant  depth  plate  girder  and  specifies  all 
of  the  details  of  the  design,  including  the  flange  thick- 
nesses, the  flange  splice  locations  and  the  locations  of  the 
transverse  and  longitudinal  stiffeners  if  they  are  required. 

1  .  5    Scope  qj_   the  I  nves  t  i  gat  ion 
Since  the  GAD- I  program  has  been  developed,  it  has  been 
possible  to  design  the  optimum  welded  plate  girder  for  a 
highway  bridge.   This  girder,  however,  is  not  necessarily 
the  actual  optimum  design  for  the  problem.   There  is  a  possi- 
bility that  cither  a  composite  or  noncomposite  rolled-section 
girder  may  be  more  economical  than  the  plate  girder.   This 
is  particularly  true  if  the  spans  of  the  girder  arc  rela- 
tively short . 


in  order  to  determine  the  true  optimum  design,  a  com- 
puter-aided design  system  for  the  rolled  section  bridge 

girder  has  been  developed.   With  this  program,  a  given 
problem  can  be  designed  with  a  plate  girder  and  a  rolled- 
section  girder  and  the  true  optimum  can  be  chosen  directly. 

In  the  present  study,  the  American  Association  of   late 
Highway  Officials  bridge  code  [20]  is  used  to  control  the 
design  of  a  highway  bri'dge  girder.   The  cross  section  of  the 
girder  is  limited  to  the  36-in.  and  33-in.  wide- flange 
rolled  sections.   The  design  can  be  either  composite  or 
noncomposite  and  can  have  cover  plates  If  they  are  required. 
The  problem  is  limited  to  a  maximum  of  four  continuous  spans 
designed  for  all  types  of  static  and  Fatigue  highway  Loads. 
The  fatigue  conditions  arc  based  on  500,000  cycles  of  stress 
If  the  design  is  composite,  then  7/8- in.  diameter  shear 
connector  studs  are  provided  at  the  proper  spacing.   The 
method  of  exhaustive  search  is  used  in  the  program  because 
of  the  limited  number  of  rolled  sections  of  proper  size. 
Each  of  the  eighteen  possibilities  is  designed,  using,  the 
optimum-seeking  methods  of  Interval  halving  and  dynamic 
programming  to  determine  the  required  cover  plates. 


CHAPTER  II 

G  [RDER  ANALYSIS 


2 . 1    Pes  ign  Loads 

The  first  problem  that  confronts  the  designer  of  a 
particular  structure  is  to  determine  the  required  design- 
load  conditions  for  that  structure.   These  conditions  usual- 
ly vary  with  the  size,  purpose  and  location  of  the  structure 
The  actual  conditions  are  often  designated  by  design  speci- 
fications and  can  be  found  in  design  codes. 

The  design  of  a  highway  bridge  is  controlled  by  the 
American  Association  of  State  Highway  Officials  code  [20]. 
This  code,  known  as  the  AASIIO  code,  states  that  all  of  the 
girders  in  a  bridge  structure  must  be  designed  for  the 
largest  truck  load  and  the  corresponding  equivalent  lane 
load  which  the  bridge  must  sustain.   The  code  designates 
the  loads  by  the  size  of  the  truck  being  considered,  i.e., 
IIS20.   Once  the  required  designation  is  chosen,  the  actual 
loading  parameters  can  be  determined.   These  parameters 
include  the  axle  loads,  the  axle  spacings,  the  equivalent 
lane  loads  and  t  h  e  corresponding  concentrated  loads. 
Table  2.1  contains  a  summary  of  the  load  parameters  for 
all  of  the  AASHO  truck  designations. 


Table  2.1:   AASHO  Load  Designations  and  Parameters 


Truck  Loadings 

i  quiva 1 cnt  Lane 
Load  i  ng 

Axle  Loads 

Axle  Spacing s 

llni  form 

Concentrated 

Des. 

1 

7 

5 

1  -  2 

Min 
2  -  3 

Max 
2  -  3 

Moment 

Shear 

H10 

1115 
H20 

IIS15 
US  2  0 

4k 
6k 
8k 
6k 
8k 

16k 

24k 
32k 
24  k 

3  2k 

24k 

32k 

14' 
14' 
14  ' 
14  ' 
14' 

14' 
14  ' 

30' 
30' 
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If  the  bridge  is  being  designed  as  part  of  the  Inter- 
state Highway  System,  the  girders  must  also  be  designed  for 
a  specified  interstate  loading  condition.   This  loading 
condition  is  based  on  the  dimensions  and  weight  of  the 
largest  military  vehicle  which  would  use  the  structure 
during  a  national  emergency.   The  loading  consists  of  two, 
24 -kip  concentrated  loads  with  an  axle  spacing  of  4  ft. 
The  effect  of  this  loading  must  he  calculated  in  order  to 
determine  if  it  g  o v  e r  n s  t  h  e  g  i  r d e r  d e  s  i  g  n . 

Another  possible  girder  loading  condition  is  produced 
by  pedestrian  pathways.   If  the  bridge  is  supplied  with 
sidewalks,  then  the  sidewalk  loads  must  be  calculated. 
According  to  the  AASHO  code,  the  following  set  of  formulas 
is  used  in  calculating  the  sidewalk  load: 


p  =  8  5 


for    Ls    <    25    ft 


P    =    60 


for    25    ft    <    Ls    <    100    ft 


P    = 


30 


3000 


'S 


5  5  -  IV 
50 


for    Ls    >    100    ft 
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where    P  =  the  load  in  pounds  per  square  foot. 

Ls  =  the  loaded  length  of  sidewalk  in  feet. 
IV  =  the  width  of  sidewalk  in  feet. 
The  total  sidewalk  load  is  distributed  over  all  of    the 
girders  in  the  structure.   Although  the  moments  due  to  this 
loading  condition  are  usually  small,  they  are  included  to 
provide  a  complete  analysis. 

Besides  the  above  loads,  which  are  considered  to  be  the 
live  loads  for  the  structure,  there  are  also  dead  loads  which 
must  be  considered.   Dead  loads  include  the  weight  of  the 
steel  girder,  the  weight  of  the  deck  slab  and  the  weight  of 
the  bridge  railing,  light  fixtures,  and  similar  loads  placed 
on  the  constructed  bridge.   If  the  design  of  the  girder  is 
noncompos i te ,  then  all  of  these  loads  are  combined  in  one 
dead  load  condition.   If  the  design  is  composite,  however, 
the  loads  are  placed  into  two  groups  because  the}-  act  on 
different  composite  sections.   The  first  group  is  considered 
the  dead  load  and  includes  the  weight  of  the  slab  and  the 
steel  girder.   The  second  group  is  the  superimposed  dead 
load  (sometimes  referred  to  as  the  long-term  live  load). 
This  group  contains  dead  loads  applied  to  the  composite 
girder  after  the  concrete  has  reached  its  maxii  un  strength 
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capacity.   An  example  of  a  superimposed  dead  load  is  the 
weight  of  the  bridge  railings  and  the  Light  fixtures. 

2  .  2   Ana  Lysis  Theo ry 


Once  the  loads  which  act  on  the  structure  have  been 
defined,  the  stresses  caused  by  these  loads  must  be  deter- 
mined.  The  design  loads  must  be  placed  on  the  structure 
so  that  maximum  design  conditions  are  developed  at  the 
critical  points  in  the  structure.   The  determination  of 
these  design  values  constitutes  a  complete  analysis  of  the 
structure  being  designed. 

Since  the  loading  conditions,  as  shown  above,  for  the 
highway  bridge  girder  are  complex  in  nature,  the  method  of 
analysis  which  is  most  efficient  is  the  method  of  influence 
lines.   This  method  enables  the  structure  to  be  loaded  with 
a  large  number  of  loading  conditions  and  the  controlling 
maximum  to  be  chosen  from  these  conditions. 

An  influence  line  is  a  graph  showing  the  variation  in 
a  particular  function  at  a  point  with  respect  to  the  posi- 
tion of  a  unit  load.   flic  ordinates  of  the  influence  line 
arc  referred  to  as  the  influence  coefficients.   The  ordinate 
at  point  j  for  the  function  I'  at  point  i,  usually  called  the 

influence  coefficient  a. .,  is  defined  as  the  value  F.  due  to 

lj  i 

a  unit  load  at  point  j  [21],   With  this  definition  various 

theorems  can  be  developed  pertaining  to  the  use  of  influence 

lines.   Norris  and  Wilbur  [--]  describe  the  following  four 
basic  Influence  line  theorems: 
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Theorem  1  -  The  maximum  value  of  a  function  due  to 
a  single  concentrated  load  is  found  by 
placing  the  load  over  the  maximum 
influence  coordinate. 

Theorem  2  -  The  value  of  a  function  due  to  a  single 
concentrated  load  equals  the  product  of 
the  load  and  the  influence  coefficient 
at  the  load  location. 

Theorem  2  is  stated  mathematically  as 

F .  =  P .  x  a . .  2.2 

Prom  Theorem  1,  the  maximum  value  of  P.,  F-;    ,  is  found 
'  l    ■'■max ' 

when  the  single  concentrated  load  is  placed  at  the  point 
where  a. .  is  the  maximum.   By  using  superposition,  the  value 

1J  .        t.     r    t- 

of  the  function  F  due  to  a  series  of  concentrated  loads  can 
be  found  from  Equation  2.5. 


F  .  =  J      P  .  x  a  .  . 


m 

I 
j  =  l 

where   m  is  the  number  of  loads. 

Theorem  3  -  The  maximum  value  due  to  a  distributed 

load  occurs  when  the  structure  is  loaded 
over  those  portions  which  have  influence 
coefficients  with  the  same  sign  as  the 
character  ol~    the  function  desired. 

Theorem  4  -  The  value  of  the  function  due  to  a  dis- 
tributed load  is  the  product  of  the  load 
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magnitude  and  the  net  area  under  the 
portion  of  the  influence  line  loaded. 
Theorem  4  is  stated  mathematically  as 

g(x)  •  a(x)dx 


2.4 

where   r  and-s  are  the  load  limits. 
The  derivation  of  Equation  2.4  can  he  found  in  the  structural 
analysis  texts  by  Norris  and  Wilbur  [22]  and  by  Shedd  and 
Vawter  [ 2  3 ] . 

There  are  two  simplifications  that  can  be  made  to  liqua- 
tion 2.4  when  dealing  with  highway  structures.   In  all  cases, 
the  distributed  load,  g(x),  can  be  replaced  by  a  constant 
representing  the  uniformly  distributed  loading  condition. 
The  uniform  load  is  usually  represented  by  w.   The  second 
simplification  deals  with  the  influence  line  function,  a(x). 
Since  this  function  is  usually  fairly  complex  in  nature  and 
too  difficult  to  obtain,  the  coefficients  are  determined  at 
a  number  of  discrete  points  and  the  values  in  between  are 
found  using  linear  interpolation.   The  value  of  the  function 
F,  due  to  a  uniform  loading  condition,  over  a  pieccwise 
continuous  influence  line,  can  be  determined  by  using  Equa- 
tion 2.5. 


Fi  =  J      llai,k  +  ai,k+1jH 

K=  1 


2.5 


where   n  is  the  number  of  elements 
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With  the  theorems  and  the  equations  expressed  above, 
any  effect  of  a  Loading  condition  can  be  obtained,  provided 
the  influence  coefficients  can  be  calculated.   The  influence 
line  is  developed  by  computing  the  required  influence  coeffi 
cients  for  each  position  of  a  unit  concentrated  load.   The 
process  is  quite  simple  for  statically  determinate  girders', 
but  becomes  much  more  complicated  for  indeterminate  or  con- 
t  inuous  g  irders . 

For  the  statically  determinate  problem,  all  that  is 
required  to  solve  for  any  influence  coefficient  are  the 
equilibrium  equations.   The  unit  load  is  placed  on  the 
structure  and  the  reactions  are  determined  using  the  equa- 
tions shown  as  liquation  2.6: 

Fm,   r+      ,  =  o 

-  left  end 

2.6 

L   vertical 

Once  the  reactions  are  found,  the  influence  coefficient  for 
any  function  can  be  determined. 

If  the  girder  being  designed  is  statically  indetermi- 
nate, then  equations  of  compatibility  arc  required.   One 
compatibility  equation  is  needed  for  each  redundant  reaction. 
These  equations  refer  to  the  requirement  that  all  the  dis- 
placements throughout  the  structure  are  consistent.   Killems 
and  Lucas  [211  note  that  the  compatibility  conditions 
require  that  the  displacements  at  supports  must  be  externally 
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satisfied.   In  the  particular  case  of  unyielding  supports, 
the  total  displacement  must  be  equa]  to  zero. 

The  compatibility  equations  can  be  developed  with 
reference  to  figure  2.1.   The  structure  is  made  statically 
determinate  by  removing  the  interior  supports  (Figure  2.1b) 
With  the  determinate  structure  defined,  the  deflections  at 
the  interior  support  points,  due  to  the  external  loads  on 
the  simple  girder,  can  be  determined.   Next  the  flexibility 
coefficients  are  determined  (Figure  2.1c  and  d)  and  the 
equations  can  be  formed: 


R2f22  +  R3f32  -  n2  =  n 


Mn  +  R3f33  -  D3  -  ° 


2.7a 


In  matrix  form,  these  equations  may  be  written 


f     f 

1  2  2     x  3  2 


f       f 

2  3      13   J 


Ft., 


3   J 


2.  7b 


or 


f  F  ] 


[  R  ]   =  [  D 
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The  F  matrix  is  called  the  flexibility  matrix  and  is  inde- 
pendent of  the  external  loads  on  the  structure.   The  coeffi- 
cients of  this  matrix  represent  the  deflections  at  the 
interior  support  locations  due  to  unit  loads  at  the  interior 
support  locations.   These  deflections  can  be  found  by  using 
the  moment -area  theorems. 
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2.1  a:   Indeterminate  Structure 


2.1  b:   Determinate  Structure 


2.1  d :   Unit  Load  at  Point  3 


Figure  2.1:   Determination  of  Compatibility  Equations 
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The  moment-area  theorem  used  to  deteri  Lne  the  del  lec- 
tion is  called  the  second  moment-area  theorem.   The  theorem 
states  that  the  deflection  at  point  1,  relative  to  a  tangent 
at  point  2,  is  equal  to  the  moment  of  the  area  under  the 
bending  moment  diagram,  between  points  1  and  2,  about  point 
1  divided  by  the  flexural  rigidity  oi'    the  girder.   l'ippard 
and  Baker  [25]* give  the  following  mathematical  form  of  the 
theorem: 


D 


EI 


Mxdx 
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If  the  moment  of  inertia  varies  along  the  girder  and  the 
moment  diagram  is  given  as  a  piecewisc  continuous  curve, 
varying  linearly  between  the  analysis  points,  then  Equation 
2.9  takes  the  following  form: 

1_   y   1/        ^  l 


D 


E  i-i  2 


M.+M.       ■  x. 

i   i+  i  I  I  .   i 


2.10 


where   x.  is  the  distance  from  the  center  of  the 
i 

element  to  the  center  of  moments. 

Once  the  required  deflections  are  determined  and  the 
compatibility  equations,  2. "a,  are  defined,  the  values  of 
the  redundant  reactions  can  be  determined.   With  these 
reactions  and  the  equilibrium  equations,  2.0,    the  remaining 
reactions  can  he  calculated. 

In  order  to  develop  the  influence  line,  the  unit  load 
must  he  located  at  each  point  along  the  structure.   This 
means  that  there  would  he  NA  different  matrix  equations  oC 
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the  type  shown  in  Equation  2.S,  where  NA  is  the  number  of 
analysis  points.   Since  the  flexibility  matrix,  however, 
docs  not  change  with  the  location  of  the  external  load,  the 
only  part  of  the  equation  that  wil]  vary  is  the  D  matrix. 
For  this  reason  it  is  possible  to  set  up  one  matrix  of  NA 
columns  for  all  the  positions  of  the  unit  load.   Also  by 
adding  two  more  rows,  which  represent  the  equilibrium  equa- 
tions, all  of  the  reaction  influence  coefficients  can  be 
determined  at  the  same  time.   Shown  in  Figure  2.2  is  a 
three  -span  girder  with  seven  analysis  points.   The  matrix 
equation  required  to  solve  for  the  reaction  influence 
coefficients  is 


[  CC  ]   [  R  ] 


where 


[  CC  ] 


[  RIL  ] 


11  1  1 

0   Lj  L2  L3 

0    d  3  3  d  5  3  0 

0   d35  d55  0 


.11 


an. 


[  Rib  ] 


1 

1 

1 

1 

1 

1 

1 

Xi 

x2 

x3 

X  ,, 

X5 

x6 

x7 

d31 

d32 

d33 

d3. 

d35 

d36 

d3 

5  1 


5  2 


5  3 


5  5 


5  6 


3  7 


5  7-1 


This  equation  is  the  same  as  the  one  which  appears  in  the 
analysis  development  by  Goble  and  DeSantis  [18]. 

The  solution  to  this  matrix  equation  can  be  approached 
in  a  number  of  ways.   The  most  efficient  method  is  one  in 
which  the  final  solution  for  the  reaction  influence  lines 
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Figure  2.2:   Sample  Girder  for  Equation  .Pevel  opment 


is  found  in  the  RIL  matrix.   By  using  a  Gaussian  search 
technique,  the  CC  matrix  is  reduced  to  an  identity  matrix 
by  simple  row  and  column  operations  common  to  matrix 
algebra.   By  repeating  the  same  operations  on  the  RIL 
matrix,  the  reaction  influence  coefficients  are  determined 
and  stored  in  the  RIL  matrix.   A  complete  description  of 
the  Gaussian  search  technique  employed  in  the  solution  is 
given  by  Coble  and  DeSantis  [19]. 

With  the  solution  for  the  reaction  influence  coeffici- 
ent complete,  all  of  the  remaining  influence  lines  for  the 
complete  analysis  can  be  calculated  by  simple  statics. 
The  influence  lines  can  then  be  loaded  and  all  the  design 
conditions  for  the  problem  can  be  determined.   With  these 
conditions,  the  girder  can  be  designed  and  optimized. 

All  of  the  theorems  and  equations  discussed  above  arc 
used  to  determine  the  design  reactions,  moments  and  shears 
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After  a  fiiKil  design  has  been  chosen,  the  \ASHO  code 
requires  that  the  deflections  be  calculated  and  compared 
to  certain  allowable  live-load  deflections.   In  order  to 
calculate  the  maximum  deflections,  the  deflection  influence 
coefficients  must  be  determined.   The  amount  of  calcula- 
tions, however,  can  be  reduced  by  employing  Maxwell's  law 
of  reciprocal  deflections.   The  law  says  that,  on  a  girder, 
the  deflection  at  point  1  due  to  a  unit  load  at  point  2  is 
equal  to  the  deflection  at  point  2  due  to  a  unit  load  at 
point  1.   This  enables  the  deflection  influence  coefficients 
for  a  point  j  to  be  determined  by  calculating  all  of  the 
deflections  along  the  girder  due  to  a  unit  load  at  point  j. 
The  maximum  deflections  can  then  be  calculated  by  using  the 
influence  coefficients  and  the  design  loads.   These  maximums 
can  then  be  compared  to  the  allowable  deflections  specified 
by  the  AASHO  code. 


2 . 3    Analysis  Program 
The  analysis  portion  of  the  computer  program  uses  the 
theory  described  in  Section  2.2  to  develop  the  values 
required  by  the  design  portion  of  the  program.   A  basic 
flowchart  for  the  analysis  program  is  shown  in  Figure  2.5. 
The  step  numbers, used  in  the  following  program  description, 
refer  to  the  numbers  in  parentheses  on  the  flowchart. 


ViD 

Enter  the  ANAL 
Subroutine 

I 


(2) 


Enter  the  REAC 
Subroutine 


I  =  1 


I  >  NA 


I  =  1  +  1 


True 


-M  a 


F  a  1  s  e 


Determine  the  simple 
span  bending  moments 


I 


(3) 


Enter  the  EQSFT  subroutine  and 
determine  the  simple  span  deflections 


I 


(4) 


Set  up  the  cut i  re 
K 1 L  matrix 


Figure  2.5:   Analysis  Program  Flowchart 
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Set  up  the 
CC  matrix 


I  {6} 


Use  the  SOLVE  routine  to  determine 
the  reaction  influence  coefficients 


Load  the  reaction 
influence  line 


I 


(9) 


Determine  the 
design  reactions 


Xo 


Yes 


-© 


Fi  eure  2.3:   ("Cont  '  d  . 


© 


{(13  ) 

Determine  the  moment  influence 
lines  and  load  the  lines 


I 


Determine  the 
design  moments 


No 


Determine  and  load  the 
shear  influence  lines 


I 


Determine  the 
design  shears 


No 


Yes 


Figure  2.5:   (Cont'd.) 
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Step  1 
The  analysis  program  is  started  by  calling  the  ANAL 
subroutine.   This  subroutine  calls  on  the  other  subroutines 
in  the  proper  order  required  to  produce  a  set  of  usable 
design  conditions. 

Step  2 
The  REAC  subroutine  is  called  on  to  solve  for  the 
reaction  influence  lines.   The  simple-span  girder,  defined 
as  the  problem  girder  with  the  interior  supports  removed, 
is  loaded  with  a  unit  load  at  each  analysis  point,  and  the 
bending  moments  are  determined.   These  simple-span  bending 
moments  are  transferred  to  the  EQSET  subroutine. 

Step  3 
The  EQSET  subroutine  is  used  to  calculate  the  simple- 
span  deflections  at  the  interior  support  points.   The  value 
of  the  tangent  at  the  left  end  (TALE)  is  determined  by  using 
Equation  2.10  and  is  then  used  in  the  following  equation  to 
determine  the  required  deflections: 


d.  =  TALI',  x  COOR(I)  -  WA(I)  2.12 


Step  4 
The  RIL  matrix  is  prepared  by  using  the  deflections 
obtained  in  step  3  and  the  equilibrium  equations.   The  RIL 
matrix  is  the  right-hand-side  matrix  of  liquation  2.11. 


Step  5 
The  coefficient  matrix  (CC)  is  prepared  by  using  the 
appropriate  values  of  the  deflections  found  in  step  3  and 
the  coefficients  of  the  equilibrium  equations.   This  CC 
matrix  corresponds  to  the  matrix  CC  in  Equation  2.11. 

Step  6 
The  SOLVE  subroutine  is  called  to  determine  the  solu- 
tion to  Equation  2.11.   By  using  a  Gaussian  search  technique, 
the  CC  matrix  is  reduced  to  an  identity  matrix  and  the  reac- 
tion influence  coefficients  are  determined  and  stored  in 
the  RIL  matrix. 

Step  7 
Steps  2  thru  6  are  repeated  using  a  different  moment 
of  inertia  for  each  subelement.   By  doing  this,  the  differ- 
ent influence  coefficients  for  the  steel  section,  the  low 
modulus  concrete  section,  and  the  high  modulus  concrete 
section  are  determined.   Bach  set  of  coefficients  is  stored 
in  its  proper  location  in  the  RIL  matrix.   Control  is  now 
passed  back  to  the  ANAL  subroutine  where  the  reaction 
influence  coefficients  are  printed  if  IPT10  is  greater 
than  2. 


Step  8 
The  reaction  influence  line  is  Loaded  by  calling  on 
the  LOAD  subroutine  and  the  values  of  the  ten  [10]  loading 


conditions  are  returned  from  this  subroutine.   Since  the 
LOAD  subroutine  is  used  extensively  in  the  analysis,  it 

is  described  in  Section  2.4. 


Step  9 
The  design  reactions  (DESREA)  are  now  determined.   The 
dead-load  reaction  is  set  equal  to  ST0REC1)  and  the  super-  ' 
imposed  dead-load  reaction  is  set  equal  to  STORE(IO).   The 
positive  and  negative  live- load  design  reactions  are  equal 
to  the  maximum  live- load  reactions  possible.   These  maximums 
are  obtained  by  combining  the  sidewalk  live  load  with  the 
other  individual  live  loads  of  the  same  sign  until  the 
largest  value  of  the  sums  is  found.   If  the  sidewalk  live 
load  is  zero,  then  the  positive  and  negative  design  reactions 
are  respectively  equal  to  the  largest  positive  and  negative 
truck  load,  lane  load  or  interstate  load. 

Step  10 
Steps  8  and  9    arc  repeated  for  each  support  along  the 
girder.   The  maximum  number  of  repetitions  is  NS.   The 
DESREA  array  is  now  printed  out  if  the  1PT10  parameter  is 
greater  tha n  zero . 

Step  1 1 
The  reaction  influence  coefficients  are  used  to  deter- 
mine the  moment  influence  line  at  a  particular  point.   This 
influence  Line  is  then  loaded  by  the  LOAD  subroutine  and  the 
ten  loading  conditions  arc  determined.   The  design  moments 


moments  (DESMOM)  are  determined  in  the  same  manner  as  the" 
design  reactions  in  step  9. 

Step  12 
Step  11  is  repeated  for  each  analysis  point  along  the 
structure  and  the  DESMOM  array  is  printed  if  N'TIO  is 
greater  than  0.   The  maximum  number  of  repetitions  is  NA . 

Step  15 
The  reaction  influence  coefficients  are  used  to  deter- 
mine the  shear  influence  line  and  this  line  is  loaded  by 
the  LOAD  subroutine  to  calculate  the  values  of  the  ten 
loading  conditions.   These  conditions  are  used  to  calculate 
the  design  shears  (DESSI1)  in  the  same  manner  as  described 
in  step  9. 

Step  14 
Step  13  is  repeated  for  the  locations  immediately  to 
the  left  and  to  the  right  of  each  analysis  point.   The  maxi- 
mum number  of  repetitions  is  2NA-2.   (The  minus  2  is  neces- 
sary because  there  is  only  one  repetition  made  for  the 
first  and  the  last  analysis  points.)   The  DESSI1  array  is 
printed  if  tPTIO  is  greater  than  0. 

Step  15 
The  analysis  of  the  girder  is  complete  and  the  control 
is  returned  to  the  main  program  for  design. 


2 . 4    Load  Subrouti nc 
The  load  subroutine  is  used  to  load  any  given  influence 
line  with  the  complex  loadings  required  in  a  highway  girder 
analysis.   It  should  be  noted  that  a  given  influence  line 
actually  consists  of  three  different  influence  lines,  if 
the  design  is  composite.   There  is  one  line  for  the  steel 
section,  one  for  the  low  modulus  concrete  section,  and  one 
for  the  high  modulus  concrete  section.   All  three  lines  are 
calculated  in  the  ANAL  subroutine  and  are  transferred  to  the 
LOAD  subroutine  at  the  same  time. 

Step  1 
The  steel  section  influence  line  is  set  into  the  RO 
array  and  the  properties  of  this  line  are  determined  by  the 
ILPROP  subroutine.   The  properties  which  are  required  by 
the  LOAD  routine  are  the  positive  and  negative  areas  and 
the  locations  of  the  maximum  and  minimum  ordinates. 

Step  2 
The  unit  dead  load  is  calculated.   Using  the  area 
properties  of  the  influence  line,  the  dead- load  loading  is 
found  by  using  Equation  2.5  with  the  load  extending  over 
the  entire  girder.   The  dead -load  loading  is  then  stored 
in  STORE! 1 ) . 
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If  the  girder  is  noncompositc,  then  the  .superimposed 
dead-load  loading  (STORE(IO))  is  set  equal  to  zero.   If  the 
design  is  composite,  then  the  low  modulus  concrete  section 
influence  line  is  set  into  the  RO  array  and  the  properties 
are  found.   The  superimposed  dead  load  is  determined  and 
the  values  of  STORE(10)-is  calculated  using%  Equation  2.5 
with  the  load  extending  over  the  entire  girder. 

Step  4 
If  the  design  is  noncompositc,  then  the  live-load 
determination  can  be  started  immediately.   If  the  design 
is  composite,  however,  the  high  modulus  concrete  section 
influence  line  is  set  into  the  RO  matrix  and  the  properties 
are  determined  by  the  ILPROP  subroutine. 

Step  5 
The  IMPACT  subroutine  is  called  on  to  determine  the 
impact  factor  required  by  the  AASHO  code.   This  factor  (HT) 
is  determined  from  the  following  formula: 

HT  =  1.0  + ,S01  „r  2.13 

D I NG - 1 2 5 

The  code  places  an  upper  limit  of  1.5  on  the  impact  factor. 

Step  6 
The  positive  lane  load  is  found  by  loading  the  positive 
area  with  the  equivalent  lane  load  and  by  placing  the 
required  concentrated  load  over  the  maximum  ordinate.   The 


load  is  Increased  by  the  Impact  factor  and  stored  in 
ST0RE(2).   The  negative  lane  load  is  obtained  by  loading 
the  negative  area  with  the  equivalent  lane  load  and  by 
placing  the  required  concentrated  load  over  the  minimum 
ordinate.   The  AASHO  code  also  requires  a  second  concen- 
trated load  placed  in  one  other  span  as  to  produce  the 
maximum  negative  moment.   The  negative  lane  load  is  then 
increased  by  the  impact  factor  and  stored  in  ST0RE(3). 

Step  7 
The  sidewalk  live  load  is  found  by  multiplying  the 
uniform  sidewalk  load  by  the  positive  and  negative  areas 
respectively.   The  positive  sidewalk  load  is  placed  in 
STORE  (8)  and  the  negative  sidewalk  load  is  placed  in  STORE  (9  J 

Step  8 
If  the  interstate  loading  is  desired,  the  influence 
line  is  loaded  with  the  two  interstate  concentrated  loads 
which  are  spaced  4  ft  apart.   There  arc  three  loading  con- 
ditions considered  in  the  interstate  loading.   The  first 
two  have  one  load  on  the  maximum  ordinate  and  the  second 
4  ft  to  the  right  or  4  ft  to  the  left.   The  third  condition 
has  the  two  loads  centered  about  the  maximum  ordinate.   The 
total  maximum  condition  is  increased  by  the  if. pact  factor 
and  stored  in  ST0RE(4).   To  find  the  negative  interstate 
load,  the  sane  three  conditions  are  repeated  using  the 


minimum  ordinate  as  the  reference.   The  negative  interstate 
load  is  stored  in  ST0RE(5). 

Step  9 
The  truck  load  analysis  begins  by  moving  the  truck 
from  left  to  right.   A  complete  analysis  consists  of  placing 
each  of  the  truck's  axles  over  the  point  in  question  and 
calculating  the  loads  due  to  the  combinations  of  the 
remaining  axles. 

Step  10 
The  first  step  is  to  place  the  rear  axle  on  the 
analysis  point  and  to  use  the  maximum  spacing  for  the  front 
axles.   The  values  of  the  influence  line  are  determined  by 
using  the  interpolation  subroutine  (ILINT).   The  resultant 
load  is  found  by  using  Equation  2.5.   This  load  is  compared 
to  the  current  maximums.   (The  positive  maximum  is  found 
in  POS  and  the  negative  maximum  is  found  in  SEG . )   The 
calculation  is  repeated  using  the  minimum  spac.ings  for  the 
front  axles. 

Step  11 
Next  the  front  axle  is  placed  on  the  point  in  question 
and  the  second  axle  is  located  the  specified  fixed  distance 
away.   The  effect  of  these  two  axles  is  calculated  and  com- 
pared with  40  percent  of  the  current  maximums.   It  is 
assumed  that  if  the  effect  of  the  front  axles  is  less  than 
40  percent  of  the  maximums,  then  the  effect  o1'    the  rear 
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axle  can  be  neglected.   If  the  effect  of  the  rear  axle  must 
be  calculated,  the  rear  axle  is  placed  at  the  minimum 
distance,  at  the  maximum  distance,  and  at  each  analysis 
point  in  between.   Goble  and  DeSantis  [18  and  19]  refer  to 
this  as  "wiggling  the  rear  axles."   Each  total  effect  is 
compared  to  the  current  maximums  and  stored  in  POS  and  SEG 
if  required. 

Step  12 
The  final  loading  for  a  particular  point  is  to  place 
the  second  axle  on  the  point  and  the  first  axle  the  speci- 
fied fixed  distance  ahead.   The  effect  of  these  two  loads 
is  calculated  and  compared  in  the  same  manner  as  in  step  11. 
If  necessary,  the  rear  axle  is  wiggled  and  the  total  effects 
are  compared  to  the  current  maximums . 

Step  15 
The  truck  load  is  repeated  for  the  truck  moving  from 
right  to  left.   The  final  maximums  are  increased  by  the 
impact  factor  and  the  positive  load  is  stored  in  ST0RE(6) 
and  the  negative  load  is'  stored  in  STORE  [7). 

Step  14 
The  STORE  array  is  now  returned  to  the  ANAL  subroutine 
where  it  is  used  to  determine  the  design  loading  conditions. 
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CHAPTER  J  1  I 
GIRDER  DESIGN 

3.1    Girder  Design  Theorv 
With  the  design  conditions  defined  at  each  analysis 
point,  the  cross  section  can  now  be  proportioned  to  carry 
the  necessary  design  stresses.   If  the  design  is  noncompos i te , 
then  the  simple,  pure  bending  theory  can  be  used  to  determine 
the  stresses.   The  equation  which  expresses  this  pure  bending 
theory  is  known  as  the  flexure  formula  and  is  shown  as  Equa- 
tion 3.1: 

f  =  M.  .  x  -  3.1 

tot    i 

or 

f  =  Mtot/S 
where   S  =  the  section  modulus  of  the  section  and 
is  equal  to  I/c . 
A  derivation  of  this  formula  is  presented  by  Byars  and 
Snyder  [26]  and  many  other  mechanics  textbooks. 

The  section  used  in  a  noncompos i te  design  is  symmetrical 
about  the  horizontal  axis.   This  al  loves  the  top  and  bottom 
fiber  stresses  of  the  steel  section  to  be  equal.   The  con- 
crete in  a  noncomposite  design  is  assumed  to  carry  no  longi- 
tudinal stresses  because  there  is  no  provision  for  the 
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transfer  of  shear  between  the  slab  and  t lie  beam.  If  pro- 
vision is  made  for  this  transfer,  then  the  design  becomes 
compos  i  tc . 

The  design  of  a  composite  bridge  girder  takes  advantage 
of  the  high  compressive  strength  available  in  the  concrete 
slab.   Some  of  the  advantages  of  composite  action  expressed 
by  Viest,  Fountain,  and  Singleton  [27]  and  by  McCormac  [28] 
are : 

1.  A  savings  in  steel  costs  is  realized, 

2.  The  beam  depth  can  be  reduced, 

3.  A  larger  percentage  of  steel  is  in  tension, 

4.  The  deflections  are  reduced  due  to  the  increase 
in  the  beam  stiffness, 

5.  A  composite  section  can  withstand  a  greater 
overload  than  a  noncomposite  design. 

The  one  disadvantage  of  the  composite  action  is  the  cost  of 
providing  the  shear  connectors.   This  disadvantage,  however, 
only  takes  precedence  over  the  advantages  in  short,  lightly 
loaded  structures. 

In  designing  a  composite  girder  by  the  AASHO  code,  the 
properties  of  the  steel  beam  and  slab  are  established  on  the 
basis  of  the  moment  of  inertia  of  the  composite  section. 
Bresler,  Lin  and  Scalzi  [29]  give  the  following  two  basic 
assumptions  used  in  the  analysis  and  design  of  composite 
beam s  £ o r  h i  g hw ay  g i r d e r s : 

1.   The  slab  is  connected  throughout  the  length 
o  f  tli  e  g  i  r  d  e  r  , 
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2.   The  stress  and  the  strain  are  linear  across 
the  depth  of  the  member. 
The  AASHO  code,  unlike  the  American  Institute  of  Steel 
Construction  specification  [30],  docs  not  allow  any  slippage 
for  the  shear  connectors. 

The  composite  design  section  shown  in  Figure  3.1  is 
composed  of  a  rolled  section,  an  effective  slab  section, 
a  bottom  cover  plate,  and  a  haunch.   The  section  could  also 
be  designed  without  the  bottom  plate.   It  is  also  possible 
for  the  slab  to  be  on  one  side  of  the  girder.   This  type  of 
section  would  be  used  for  a  facia  girder  and  is  shown  in 
Figure  3.2. 

The  effective  slab  section  is  defined  as  a  section  of 
the  bridge  deck  having  an  effective  slab  width  and  a  thick- 
ness equal  to  the  structural  depth  of  the  slab.   The  AASHO 
code  defines  the  effective  slab  width  of  the  section  in 
Figure  3.1  as  not  exceeding  any  of  the  following: 

1.  One- fourth  of  the  girder  span, 

2.  The  distance  center- to- center  of  girders, 

3.  Twelve  times  the  slab  thickness. 

If  the  composite  section  is  similar  to  the  facia  girder 
shown  in  Figure  3.2,  then  the  effective  slab  width  should 
not  exceed  any  of  the  following: 

1.  One- twelfth  of  the  span  length, 

2.  One-half  the  center- to-center  distance  to 
the  adj acent  gi rder , 

3.  Six  times  the  slab  thickness. 
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Figure  3.1:   Standard  Composite  Girder  Cross  Section 
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Figure  5.2:   Facia  Girder  Composite  Cross  Section 
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The  section  properties  of  the  composite  section  are 
determined  by  the  transformed  area  method.   This  method 
assumes  that  the  bond  between  the  concrete  and  steel  is 
strong  enough  so  that  the  strains  in  the  concrete  and  steel 
at  the  junction  are  equal.   Due  to  this  assumption,  the 
relationship  between  the  stresses  at  a  given  distance  from 
the  neutral  axis  is 


where 


-  1 

f   =  -  f 

c    n   s 


n  =  E  /E 
s   c 


3.2 


The  factor  n  is  called  the  modular  ratio.   The  area  of  the 
concrete  slab  is  transformed  into  an  equivalent  area  of 
steel  by  dividing  the  effective  slab  width  by  the  modular 
ratio.   The  values  of  the  modular  ratio  specified  by  the 
AASHO  code  are  shown  in  Table  3.1. 

Table  5.1:   Values  of  the  Modular  Ratio 


Concrete 
Strong  tli 

n 

2000  -  2400 
2500  -  2900 
3000  -  5900 
40  0  0  -  4900 
5  0  00  or  more 

15 

12 

10 

8 

6 
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The  actual  form  of  the  flexure  formula  used  in  a  com 

positc  design  contains-  three  terms  instead  of  just  one  as 

i n  Equat ion  3.1: 

Mn.  y    M   .  y,      M    \- 
,-     DL  's     SLL  'hnc     LL  -  hmc 


I 


st 


I 


lmc 


T 


3.3 


hmc 


The  first  term  in  Equation  3.3  represents  the  stress  at  a 
point  due  to  the  dead  load  of  the  structure.   This  dead  load 
is  carried  by  the  steel  section  alone.   The  properties  of 
the  steel  section  are  independent  of  the  slab  dimensions. 

The  second  term  in  the  equation  represents  the  stress 
due  to  the  superimposed  dead  load.   This  load  must,  accord- 
ing to  the  AASHO  code,  be  carried  by  the  low  modulus  com- 
posite section.   The  properties  of  this  section  are  calcu- 
lated by  using  a  concrete  area  transformed  by  a  modular 
ratio  equal  to  5n.   The  increased  ratio  is  used  to  account 
for  creep  effects. 

The  third  term  in  the  equation  represents  the  stress 
due  to  the  live  load  on  the  structure.   The  live  load,  since 
it  is  placed  on  the  structure  after  the  concrete  has  reached 
its  maximum  strength,  acts  on  the  high  modulus  composite 
section.   The  properties  of  this  section  are  based  on  a 
concrete  area  transformed  with  a  modular  ratio  of  n. 

By  using  Equation  3.1  for  a  noncomposite  design  and 
Equation  5.3  for  a  composite  design,  the  steel  girder  cross 
section  can  he  proportioned  so  that  the  resulting  flexural 
stresses  are  within  the  allowables  required  bv  the  design 
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code.   [f  the  design  is  composite,  then  the  compressive 
stress  in  the  extreme  concrete  fiber  must  also  be  within 
the  accepted  code  allowable. 

With  the  girder  designed  for  flexural  stresses,  the 
shear  stresses  must  be  determined.   The  horizontal  and 
vertical  shear  stresses  at  any  point  in  the  cross  section 
can  be  found  from  Equation  3.4: 


v  = 


VQ 
bl 


3.4 


where   b  =  the  width  of  the  section  where  the 
shear  stress  is  desired. 
Byars  and  Snyder  [26]  derive  this  formula  and  also  demon- 
strate that,  for  an  I-shapcd  section,  the  shear  stress 
across  the  web  can  be  assumed  to  be  uniform.   Since  the 
stress  is  uniform,  the  AASHO  code  states  that  all  of  the 
external  vertical  shear  is  carried  by  the  steel  girder  web 
and  can  be  calculated  using  Equation  3.5: 


v  =  V.  ./A   , 

tot   web 


3.5 


For  this  calculation,  the  area  of  the  web  is  equal  to  the 
product  of  the  overall  depth  of  the  rolled  section  and  the 
web  t  h  ickness . 


5  .  2    Fatigue  P> e s  i  g n 

In  a  highway  bridge  girder  there  is  a  fluctuating 

stress  level  at  every  point  because  of  the  alternating 

loading  and  unloading  of  the  structure  due  to  the  live 
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loads.   This  continuous  fluctuation  causes  dislocations  in 
the  crystalline  structure  of  the  steel,  and  eventually 
causes  cracks  and  other  localized  failures.   This  type  of 
failure  is  referred  to  as  a  fatigue  failure. 

The  fatigue  failures  in  a  structure  nay  occur  at  loads 
which  are  well  below  the  allowable  static  loading  conditions, 
because  the  stress  is  continually  fluctuating  over  a  particu- 
lar range.   The  larger  the  fluctuation,  the  worse  the  condi- 
tion.  The  fatigue  condition  becomes  extremely  critical 
when  the  stress  alternates  between  tension  and  compression. 

The  number  of  fluctuations,  or  cycles,  for  which  a 
girder  can  be  designed  varies  from  100,000  cycles  to 
2,000,000  cycles.   The  exact  number  of  cycles  used  depends 
on  the  use  and  location  of  the  bridge. 

The  method  used  in  fatigue  design  is  to  reduce  the 
allowable  stresses  at  critical  points.   The  AASHO  code 
provides  the  following  formulas  for  this  fatigue  reduction: 

k,  f 
r    1  -  k?R 


and 


w  here 


.  S  5  F 
r         .55  1- 

*-  l   ro 

F 


and   R  =  the  algebraic  ratio  of  the  minimum  stress 
to  the  maximum  stress. 


The  k. ,  f   ,  and  a  are  the  fatigue  constants.   These  fatigue 

•    ro  ° 

constants  have  been  determined  primarily  from  fatigue  tests. 
A  description  of  some  of  the  tests  and  the  results  are  given 
in  the  II.  S.  Steel  bridge  manual  [51]. 

The  values  of  the  fatigue  constants  depend  on  the  num- 
ber of  cycles  of  stress,  the  category,  type  and  location  of 
material  and  the  yield  point  of  the  base  metal.   The  values 
required  for  a  particular  fatigue  calculation  can  be  found 
in  the  AASHO  code.   The  constants  for  500,000  cycles  of 
stress  are  given  in  Table  5.2.   These  are  the  values  which 
are  used  in  the  design  program. 

Table  5.2:   Fatigue  Constants 


AASHO 

Cate- 
gory 

Type  and 
Location 
of  Metal 

Type  of 
Maximum 
Stress 

Equation 

500,000  cycles 

f 
ro 

a 

k2 

A 
D 

F 

C 

Base 
Metal 

Metal 
Adj  acent 
to  Butt 
Weld 

Metal 
Adj  acent 
to  Fillet 

Weld 

Weld 
Metal 

Tension 
Compr . 

Tens  ion 
Compr . 

Tens  ion 
Compr . 

Shear 

5.6 

5.7 

5.6 

5.7 

5.  6 
5.6 

3  .  6 

20500 
13500 

17200 
106  00 

12000 
12000 

108  00 

.78 
.78 

.  2  5 
.  25 

0.0 
0.0 

.  3  6 

'.  5  5 
.62 

1  .0 
1  .0 

.55 
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hi  the  design  of  a  rolled  section,  fatigue  is  ;i  major 
factor  in  the  cover  plate  design.   The  ends  of  the  cover 
plates  are  determined  by  using  fatigue  calculations  for  the 
base  metal  adjacent  to  a  fillet  weld  [Category  1  J.   The 
allowable  stresses  in  this  case  are  greatly  reduced,  thus 
requiring  that  the  plates  be  cut  off  in  regions  of  low 
moments.   The  remaining  categories  are  checked  in  the  design, 
but  they  usually  do  not  govern  the  design  of  the  girder. 

5.3    Design  Details 


After  the  cover  plates  required  for  a  particular  rolled 
section  are  determined,  the  details  of  the  design  must  be 
calculated.   These  design  details  include  the  determination 
of  the  weld  sizes  for  the  cover  plates,  the  design  of  the 
required  bearing  stiffeners  and  the  calculation  of  the  shear 
connector  spacing  for  a  composite  girder.   In  the  optimiza- 
tion problem,  each  of  these  details  must  be  completed  for 
each  design  because  they  add  a  substantial  cost  to  the 
girder. 

The  welds  which  must  be  determined  for  each  cover  plate 
include  the  end  welds  and  the  seal  welds.   The  end  welds  arc 
designed  to  carry  the  maximum  flexural  stress  developed  at 
a  fixed  distance  from  the  ends  of  the  plate.   The  code  sets 
this  distance  as  1.5  times  the  width  of  the  cover  plate. 
The  required  stress  is  determined  using  either  Equation  3.1 
or  5.5.   V.'hcn  the  stress  in  the  plate  is  known,  Equation  5.8 
is  used  to  determine  the  force  in  the  plate: 


ILL     tot     pi 


When  the  Force  in  the  plate  is  known,  Blodgett  [31  J 
shows  that  the  size  of  the  required  fillet  weld  is 


3.8 


t  = 


Force 


707  L   £  ,, 
w   all 


3.9 


In  the  ahove  equation  the  L   is  the  total  length  of  the  weld 
1  w 

which,  according  to  the  AASIIO  code,  is  four  tines  the  plate 

width.   The  allowable  weld  metal  stress  (f  ,,)  is  the  lesser 

all 

value  of  the  allowable  static  shear  stress  and  the  allowable 
fatigue  stress,  calculated  using  Equation  3.6  and  the 
constants  listed  under  category  G  from  Table  3.2.   The 
fillet  weld  equation  then  becomes 


Force 


2.8  28  W  ,  f  . . 
pi   all 


5.10 


The  required  weld  size  is  usually  given  in  sixteenths  of  an 
inch.   The  code  requires  that  a  minimum  size  be  used  depend- 
ing on  the  thicknesses  of  the  plates  being  joined. 

The  design  of  the  seal  weld  for  a  cover  plate  is 
basically  similar  to  that  of  the  end  weld  except  that  shear 
stress  is  used  instead  of  flexural  stress.   Equation  3.4  is 
used  to  determine  the  maximum  shear  stress  developed  along 
the  cover  plate  and  then  liquation  3.11  is  used  to  determine 
the  size  of  the  seal  weld: 


1 v_ 

2  (  .  7  0  7  )  f 


al  1 


1.414  f 


5.11 


all 
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The  2  in  the  denominator  is  present  because  there  is  a  seal 
weld  on  each  side  of  the  cover  plate.   The  code  specifies 
that  the  seal  weld  must  extend  over  the  entire  length  of 
the  plate  in  order  to  prevent  corrosion  and  must  have  a 
minimum  thickness  of  5/16  in. 

The  bearing  stiffeners  are  required  to  transfer  the 
large  web  shear  stresses,  at  the  reactions,  to  prevent  web 
crippling.   In  most  of  the  wide-flange  rolled  sections, 
bearing  stiffeners  arc  not  required.   AASHO  requires  bearing 
stiffeners  in  those  girders  which  have  shear  stresses  in 
excess  of  75  percent  of  the  allowable  shear  stress.   If 
stiffeners  are  required,  they  are  designed  as  columns.   The 
column  section  used  is  composed  of  the  two  stiffener  plates 
and  a  portion  of  the  web  which  is  18  times  the  web  thickness 
in  length.   The  fillet  welds  which  connect  the  stiffener  to 
the  girder  web  are  designed  similar  to  the  seal  welds  for 
the  cover  plates. 

The  shear  connectors  for  a  composite  girder  arc  designed 
to  transfer  the  horizontal  shear  from  the  slab  to  the  steel 
section.   These  shear  connectors  are  placed  transversely 
across  the  flange  at  either  regular  or  variable  spacings. 
The  connectors  are  designed  for  fatigue  as  presented  by 
Slutter  and  Fisher  [35]  and  checked  for  ultimate  strength. 
The  spacing  at  any  location  can  be  found  as  follows: 


Spac Lng  = 


n 
r   c 
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1 1 


where   S   =  the  range  of  horizontal  shear  per  inch. 
The  allowable  shear  per  connector  depends  on  the  number  of 
fatigue  cycles  and  the  type  and  the  size  of  the  connector 
used.   The  horizontal  shear  per  inch  is  found  by  using  an 
equation  similar  to  3.4: 


S   =  V  Q/I  3.13 

r     r  K 


where   V   =  the  range  of  vertical  shear  due  to 
live  load  plus  impact. 
The  moment  of  inertia  used  in  Equation  3.13  is  the  trans- 
formed moment  of  inertia  in  the  positive  region.   If  shear 
connectors  are  placed  in  the  negative  moment  region,  then 
the  moment  of  inertia,  for  the  negative  region,  is  that  of 
the  steel  girder  and  the  longitudinal  slab  reinforcement. 

In  order  to  satisfy  the  ultimate  strength  requirements 
for  the  connectors,  a  specified  number  of  connectors  must 
be  placed  between  the  points  of  maximum  moment  and  the  points 
of  zero  dead-load  moment.   Equation  3.14  is  used  to  determine 
the  number  of  connectors  required  for  the  ultimate  strength: 


N  = 


8  5  S 
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where   I1  =  the  lesser  value  of  the  force  in  the 
steel  girder  and  the  concrete  slab. 
The  ultimate  strength  of  a  shear  connector  depends  on  the 
size  and  the  tvpe  of  connector. 
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The  final  shear  connecter  requirement  which  must  be 

satisfied  is  an  additional  number  of  connectors  to  develop 
the  slab  stress  around  the  points  of  contraf lexure.   The 

number  of  connectors  required  is 

A   f 
Nc  =  .r2  r  3.15 

r 

where   A   =  the  area  of  the  steel  reinforcement, 

and     f   =  the  range  of  live  load  stress  in  the 

reinforcement  (may  be  taken  as  10,000  psi). 

These  additional  connectors  are  needed  only  if  the  negative 

moment  region  is  noncompos i te . 

3  .  4    Design  Program 

The  design  portion  of  the  computer  program  operates  on 
each  wide-flange  rolled  section  and  completely  designs  the 
girder  for  the  given  problem.   A  complete  design  consists 
of  the  determination  of  the  cover  plates  and  their  required 
welds,  the  design  of  the  bearing  stiffeners  and  the  calcula- 
tion of  the  shear  connector  spacing. 

The  design  program  is  started  by  determining  the  lightest 
section  which  will  carry  the  loads  without  cover  plates.   This 
section  is  completely  determined  and  its  cost  becomes  the  base 
cost  for  the  optimization  problem.   The  urogram  now  enters  a 
loop  which  designs  all  of  the  remaining  rolled  sections.   A 
flowchart  for  this  loop  is  shown  in  Figure  3.3. 
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Step  1 
The  rolled  section  being  designed  is  transferred  to 
the  THDES  subroutine  where  the  thicknesses  of  the  cover 
plates  required  for  each  subelcment  are  determined. 

Step  2 
The  THDHS  subroutine  first  determines  the  maximum 
moment  in  the  subelement.  .  With  this  moment  determined, 
the  rolled  section  is  tried  without  plates.   The  individual 
stresses  are  determined  and  checked  against  the  allowable 
static  and  fatigue  base  metal  stresses.   If  these  stresses 
are  not  violated,  then  no  cover  plates  are  required  for  the 
element . 

Step  5 
If  a  cover  plate  is  required  at  a  particular  subelement, 
the  width  of  the  plate  is  determined  based  on  the  flange 
width  of  the  rolled  section.   The  thickness  of  the  plate  is 
determined  by  using  interval  halving.   The  interval  halving 
technique  uses  only  thicknesses  of  sixteenths  of  an  inch. 
The  thickness  is  varied  until  the  allowable  stress  condition 
is  satisfied.   The  thickness  and  width  of  the  required  plate 
are  stored  in  the  proper  array. 

Step  4 
Steps  2  and  3  are  repeated  for  each  of    the  subelements 
along  the  girder.   With  all  the  plate  sizes  determined,  the 


c 


;ontrol  is  passed  back  to  the  main  program 


50 


Step  5 
The  properties  of  each  subelement  are  calculated  by 
using  the  SFMIC  subroutine  and  these  properties  are  trans- 
ferred to  the  CPDES  routine. 

Step  6 
The  CPDES  subroutine  uses  the  fatigue  theory  discussed 
in  Section  3.2  to  determine  the  actual  location  of  the  ends 
of  the  required  cover  plates.   The  routine  is  used  to  deter- 
mine all  of  the  necessary  conditions  required  to  determine 
the  cutoff  point  for  each  plate.   The  actual  point  is 
located  by  the  TERMIN  subroutine. 

Step  7 
The  TERMIN  subroutine  uses  an  increasing  step  size 
optimization  technique  to  determine  the  ends  of  the  plates 
in  the  most  efficient  manner.   The  end  of  the  plate  is 
defined  as  the  point  along  the  girder  where  the  stress  docs 
not  violate  the  fatigue  allowable  stress.   >nce  this  point 
is  determined,  control  is  passed  back  to  the  CPDES  routine. 

Step  8 
Steps  6  and  7  are  repeated  until  the  fatigue  limits  oi~ 
each  cover  plate  on  the  top  and  bottom  of  the  girder  arc 
determined.   With  the  limits  found,  the  final  cover  plates 
for  the  top  and  bottom  of    the  girder  are  determined  in  the 
FLOPT  subroutine. 
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Step  9 
The  FLOPT  routine  uses  the  cost  array  and  the  cover 
plate  arrangement  to  determine  the  optimum  location  of  the 
cover  plate  splices.   A  dynamic  programming  method  is  used 
to  accomplish  this  optimization.   Upon  completion  of  this 
optimization,  the  COVPL  array  is  determined  in  its  final 
form  and  control  is  returned  to  the  main  program. 

Step  10 
The  CPWDES  subroutine  is  now  used  to  design  all  of  the 
cover  plate  welds.   The  end  welds  for  a  particular  plate  are 
designed  for  both  ends  and  the  larger  thickness  is  retained 
as  the  design.   If  the  plate  is  butt  welded  to  other  plates 
at  both  ends,  then  the  end  weld  thickness  is  set  to  zero. 
The  seal  welds  are  designed  by  the  subroutine  using  Equation 
3.11.   After  these  calculations  for  each  cover  plate  aloni* 
the  girder  are  completed,  the  control  is  returned  to  the 
main  program. 


Step  11 
The  STIFF  routine  is  called  to  design  the  bearing 
stiffeners  for  each  support.   If  the  support  does  not 
require  a  stiffener,  then  the  values  of  the  BEAR  array  are 
set  equal  to  zero.   if  a  stiffener  is  required,  the  width 
and  thickness  are  calculated  such  that  the  allowable  bearing 
and  the  allowable  compressive  stresses  arc  not  violated. 
The  welds  which  hold  the  stiffeners  in  place  are  calculated 
and  the  control  is  returned  to  the  main  program. 
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Step  12 
The  CONDES  routine  is  used  to  determine  the  spacing  of 
the  shear  connectors.   The  shear  connector  design  is  based 
on  the  use  of  7/8- in.  shear  connector  studs.   The  required 
spacing  at  each  analysis  point  is  determined  and  rounded 
off  to  the  nearest  lower  multiple  of  three  inches.   The 
spacing  is  corrected  for  the  ultimate  strength  discussed  in 
Section  3.3  and  the  CONSP  array  is  filled  with  the  required 
spacing  information.   The  routine  sets  the  least  number  of 
spaces  at  a  particular  value  at  three.   Upon  completion  of 
the  spacing  design,  the  control  is  returned  to  the  main 
program . 

Step  15 
Once  the  spacing  for  the  connectors  is  determined,  the 
design  of  a  rolled  section  is  complete.   The  OBJECT  sub- 
routine is  now  called  on  to  determine  the  weight  and,  if 
necessary,  the  cost  of  the  given  design. 

Step  14 
The  cost  or  weight  of  the  design  is  compared  to  the 
current  optimum  by  the  KEEP  routine.   If  the  girder  is  an 
improved  design,  the  elements  are  stored  in  the  proper  arrays 
If  the  girder  is  not  an  improved  design,  then  the  old  optimum 
is  retained. 
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Step  15 
Steps  1  to  ] 4  arc  repeated  for  each  of  the  remaining 
rolled  sections.   After  the  last  section  is  designed,  the 
program  takes  the  current  optimum  and  continues  into  the 
remaining  optimization  phases  of    the  program. 
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CHAPTER    [V 
PROBLEM   OPTIMIZATION 

■ 

4  . 1    Optimi zation  Theory 

Optimization,  as  stated  by  Beveridge  and  Schechter  [34], 
"is  the  collective' process  of  finding  a  set  of  conditions 
required  to  achieve  the  best  result  for  a  given  situation." 
In  other  words,  optimization  is  the  process  of  obtaining  a 
set  of  design  parameters  which  either  maximizes  or  minimizes 
a  particular  function.   The  function  being  optimized  is 
referred  to  as  the  objective  function.   This  function  may 
be  either  mathematically  simple  or  complex,  depending  on 
the  nature  of  the  problem.   The  actual  form  of  this  objective 
function  usually  determines  the  most  appropriate  method  of 
optimi  zation . 

The  design  of  a  physical  system  ideally  contains  three 
steps.   These  steps,  according  to  Wilde  and  Beightler  [35], 


are 


Determination  of  the  interaction  of  the 

system  variables, 

Development  oi'   a  simple  measure  of 

e f  fee t i vencss , 

Determination,  development,  and  solution  o 

the  most  effective  optimum-seeking  method. 
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The  first  step  is  handled  by  developing  a  set  of  constraints 
for  the  particular  problem  being  solved.   These  constraints 
can  be  imposed  by  design  codes,  tradition,  or  the  designer. 
This  step  usually  requires  the  greatest  amount  of  effort. 

The  second  step  is  the  actual  development  of  the 
objective  function.   The  formulation  of  the  function  can 
range  from  being  simple  to  being  so  complex  it  is  almost 
impossible.   The  degree  of  difficulty  in  describing  the 
function  depends  on  the  type  of  problem,  the  type  of 
effectiveness  being  measured,  and  the  degree  of  accuracy 
desired.   The  less  the  degree  of  accuracy  of  the  objective 
function,  the  less  accurate  will  be  the  developed  optimum 
solution. 

The  third  step  leads  the  designer  to  many  different 
types  of  optimization  procedures.   The  advantages  of  each 
method  vary  with  the  problem  being  solved.   It  is  the 
designer's  responsibility  to  know  the  limitations,  the 
advantages  and  the  disadvantages  of  each -method  as  they 
relate  to  the  particular  problem  being  solved. 

In  the  design  program  for  the  optimization  of  a  rolled- 
section  highway  bridge  girder,  there  are  actually  two  differ- 
ent opt  i  mum- seeking  methods  employed.   I'he  method  of 
exhaustive  search  is  used  as  the  basic  design  method  and 
dynamic  programming  is  used  to  determine  the  locations  of 
the  cover  plate  splices.   The  thicknesses  of  the  cover  plates 
are  determined  by  using  interval  halving  to  solve  the  govern- 
ing equa.t  i  on  . 
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The  exhaustive  search  technique,  sometimes  called 
exhaustive  enumeration  or  brute  force  method,  simply 
evaluates  the  objective  function  for  each  of  the  problem 
possibilities  and  picks  the  optimum  solution  directly. 
This  method  can  be  very  tedious  if  the  number  of  design 
possibilities  is  large.   It  is  therefore  used  only  when  the 
overall  design  space  for  the  problem  is  fairly  limited  in 
size. 

Since  the  wide-flange  sections  used  for  highway  girders 
include  only  the  36-  and  33-in.  sections,  there  are  only 
IS  sections  to  be  considered.   The  value  of  the  objective 
function  is  determined  for  each  wide-flange  section  after 
it  has  been  completely  designed,  including  the  cover  plates, 
the  bearing  stiffeners  and  the  shear  connectors.   The 
design  With  the  minimum  value  of  the  objective  function  is 
then  chosen  as  the  optimum  design. 

By  setting  up  the  wide-flange  sections  in  descending 
weight  order,  the  number  of  executions  in  the  exhaustive 
search  can  be  reduced.   Starting  with  the  heaviest,  and 
proceeding  in  the  descending  weight  order,  each  section  is 
tried  without  cover  plates  at  the  critical  moment  location, 
until  the  lightest  girder  which  is  acceptable  is  determined. 
By  assuming  that  the  cost  of  the  shear  connectors  and 
stiffeners  for  a  uniform  section  are  independent  of  the 
section  used,  all  of  the  sections  heavier  than  the  one  which 
has  been  determined  can  be  eliminated  from  the  search.   This 
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means  thai  a  number  of  the  rolled  sections  need  not  be 
designed  completely  and  a  resulting  savings  in  computation 
costs  i  s  real i  zed . 

Interval  halving  is  a  technique  which  can  be  used  to 
solve  equations  when  the  simple,  direct  methods  are 
inappropriate.   The  method  uses  a  bounded  interval  which 
is  continuously  reduced  in  size  until  a  certain  tolerance 
is  obtained.   By  reducing  the  interval  size  by  a  factor  of 
one-half  with  each  calculation,  the  total  number  of  neces- 
sary calculations  is  greatly  reduced  and  there  is  a  result- 
ing savings  in  calculating  costs.   The  two  necessary  condi- 
tions for  the  interval  halving  technique  are: 

1.  The  function  must  be  continuously  increasing 
or  decreasing. 

2.  The  initial  boundaries  must  be  specified. 

The  first  condition  enables  the  method  to  completely  elimi- 
nate one  half  of  the  interval  with  each  calculation,  and 
the  second  condition  enables  the  first  calculation  value  to 
be  specified. 

The  method  begins  with  the  specified  interval  and 
calculates  all  the  necessary  information  for  the  interval 
midpoint.  If  all  of  the  constraints  are  satisfied  for  this 
midpoint  value,  then  the  value  becomes  the  new  upper  bound- 
ary. This  is  possible  because  the  function  is  continuously 
increasing  or  decreasing  and  the  higher  values  of  the  vari- 
ables w i  1  1  also  satisfy  the  constraints,  but  will  cause  a 
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larger  cast.   If  the  constraints  are  not  satisfied,  then  the 
midpoint  value  becomes  the  lower  boundary.   The  new  interval 
is  used  in  the  sane  manner  as  the  initial  interval.   The 
entire  procedure  is  repeated  until  the  remaining  interval 
is  less  than  the  tolerance  desired.   Either  the  upper  or 
the  lower  bound  can  then  be  used  as  the  optimum  solution. 
In  most  problems  the  average  of  the  two  bounds  is  considered 
as  the  solution.,  provided  that  the  final  interval  is  small 
enough . 

In  determining  the  cover  plate  thicknesses  for  each 
subelement,  there  are  many  calculations  required,  including 
the  section  properties  for  each  plate  size,  the  bending 
stress  at  the  point  in  question  and  the  allowable  stresses 
for  the  fatigue  and  static  conditions.   Since,  for  each 
moment  condition  and  rolled  section  there  is  only  one  plate 
size  which  will  give  a  stress  equal  to  the  allowable  stress, 
the  function  for  the  thickness  can  be  considered  to  be  con- 
tinuously decreasing.   The  initial  boundary  conditions  can 
be  specified  by  using  the  requirements  of  the  AASHO  code. 
With  the  two  conditions  for  interval  halving  satisfied,  the' 
method  can  be  used  for  the  determination  of  the  cover  plate 
thi  cknesses . 

In  order  to  use  interval  halving,  however,  the  method 
must  be  adjusted  so  that  the  thicknesses  wi 1 1  he  determined 
in  values  of  sixteenth-inch  increments.   This  is  accomplished 
bv  adjusting  the  initial  interval  size  and  by  making  the 
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boundaries  Integers.  Each  halving  process  then  produces  an 
integer  which  corresponds  to  the  number  of  sixteenths  of  an 
inch  in  the  plate  thickness.  The  calculations  are  completed 
for  the  thickness  and  the  interval  is  reduced.  The  process 
is  continued  until  the  final  interval  size  is  1.0.  The 
optimum  plate  size  is  then  stored  in  the  proper  array  loca- 
tion. 

Dynamic  programming  is  a  method  used  in  sequential 
system  optimization  problems.   It  can  be  applied  to  situa- 
tions in  which  many  decisions  arc  required,  as  long  as  the 
decisions  made  at  later  stages  do  not  affect  the  performance 
of  earlier  stages.   A  mathematical  formulation  of  the  dynamic 
programming  method  can  be  found  in  the  optimization  texts  by 
Wilde  and  Beightler  [35],  Pierre  [36]  or  Penn  [37].   All  of 
the  mathematical  formulations  arc  based  on  Bellman's  [38] 
principle  of  optimal ity. 

Dynamic  programming  is  a  method  of  decomposition  which 
divides  the  Riven  problem  into  a  set  number  of    individual 
problems.   In  determining  the  optimum  condition  for  section 
i,  the  onlv  quantities  which  are  used  are  the  values  of  the 
variables  in  sections  i+1  and  i,  and  the  value  of  the 
objective  function  for  the  system  up  to  and  Including 
section  i-1.   The  objective  function  is  determined  for  all 
of  the  possible  combinations  of    the  variables  in  section 
i  and  i+1.   The  values  of  the  objective  ''unction  arc  then 
added  to  the  value  of  the  function  for  the  problem  through 
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section  i-1.   The  optimum  variable  can  then  be  directly 
chosen  and  placed  in  an  appropriate  array  location  defined 
by  the  section  number  i  and  the  variable  value  in  section 
i  +  1. 

The  method  begins  by  calculating  the  variable  combina- 
tions for  all  of  the  sections  and  their  adjacent  sections 
in  the  problem  sequence.   When  the  last  section  is  reached, 
the  final  values  of  the  objective  function  are  calculated. 
By  comparing  these  values,  the  optimum  for  the  last  section 
is  determined.   Then,  by  working  backwards,  all  of  the 
variables  in  each  section  can  be  determined.   This  process 
leads  to  an  optimization  of  the  entire  problem. 

The  dynamic  programming  method  is  used  to  determine 
the  final  cover  plate  combination  for  a  given  rolled  section 
The  actual  problem  is  decomposed  into  sections  corresponding 
to  the  cover  plates  which  have  been  tie  fined  by  the  fatigue 
stress  conditions.   The  variable  of  each  section  is  the 
plate  thickness,  and  the  function  being  minimized  is  the 
total  cost  of    the  cover  plate.   By  varying  the  thickness 
o  l~    each  section  and  determining  the  total  cost  of  the  plate 
by  the  procedure  described  above,  the  optimum  location  of 
the  cover  plate  splices  can  be  located.   This  method  is 
similar  to  the  flange  smoothing  method  used  by  Goble  and 
Razanj  [17].   DeSantis  and  Goble  [18]  illustrate  this 
method  for  the  optimization  ol~    the  Flanges  for  a  welded 
p  late  g  i  rder . 
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4 . 2    Ob j ect  ive  Function 
The  term  objective  function  refers  to  the  means  by 
which  the  program  evaluates  the  merit  of  a  particular  design 
For  the  rolled -sect ion  highway  girder  design,  the  function 
can  be  either  the  total  weight  or  the  total  cost  of  the 
girder.   The  designer  specifies  which  method  of  optimization 
is  desired  by  the  value  used  for  IPT6.   If  IPT6  has  a  value 
of  0,  then  the  optimization  is  based  on  cost.   The  equation 
for  the  total  girder  cost  is 


Total  Cost 


Cost  of  the  rolled  section  + 

Cost  of  the  cover  plates  + 

Cost  of  the  shear  connectors  + 

Cost  of  the  bearing  stiffeners  + 

Cost  of  the  welding. 


4.1 


If  the  value  of  IPTb  is  1,  then  the  optimization  is  based  on 
the  weight  of  the  girder.  The  equation  for  the  total  girder 
weight  is 


Total  Weight  =  Weight  of  the  rolled  section  + 

Weight  of  the  cover  plates  + 

Weight  of  the  bearing  stiffeners  + 

Weight  of  the  shear  connectors. 


4.2 


The  weight  function  first  determines  the  volume  of  steel 
used  in  each  of   the  individual  portions  of  Equation  4.2.   A 
unit  weight  of  490  pounds  per  cubic  foot  is  used  to  convert 
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the  volumes  into  weights.   These  weights  are  added  together 
to  get  the  total  weight. 

If  the  cost  of  a  girder  is  desired,  then  the  weights 
are  calculated  as  above  and  then  converted  into  costs.   The 
costs  of  the  rolled  section,  the  cover  plates,  and  the 
bearing  stiffeners  are  determined  by  multiplying  each  of 
the  weights  by  their  respective  unit  costs.   These  unit 
costs  are  expressed  in  dollars  per  pound. 

The  cost  of  the  shear  connectors  consists  of  two 
separate  elements.   The  first  element  is  the  material  cost 
of  the  connectors  and  is  based  on  the  total  weight  of  the 
connectors.   The  cost  is  calculated  by  multiplying  the  total 
connector  weight  by  the  unit  cost.   This  unit  cost  is  also 
expressed  in  dollars  per  pound.   The  second  element  is  the 
installation  cost  and  is  based  on  the  total  number  of 
connectors.   This  cost  is  calculated  by  multiplying  the 
total  number  of  connectors  by  the  installation  cost  per 
connector . 

The  welding  costs  are  determined  for  all  the  welds  on 
the  girder.   The  welds  which  are  considered  are  the  cover 
plate  butt  welds,  the  cover  plate  seal  welds,  and  the  welds 
used  to  connect  the  stiffeners  to  the  girder.   The  weld  cost 
can  be  found  for  any  weld  using  the  following  equation: 

Cweld  =  C>  +  C2(L)  +  C3^V)  4-5 
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where  C]  is  the  fixed  cost  of  the  weld,  C2(L)  is  the  vari- 
able cost  which  is  a  function  of  the  weld  length  and  C3(V) 
is  the  variable  cost  which  is  a  function  of  the  weld  volume. 

The  fixed  cost  of  a  weld  accounts  for  all  of  the  costs 
involved  in  the  set-up  time,  the  inspection,  and  the  other 
fixed  items  of  cost.   'Die  variable  cost,  based  on  the  weld 
length,  reflects  the  costs  of  x-ray  inspection  and  joint 
preparation . 

The  variable  cost,  based  on  volume,  reflects  the. actual 
cost  of  the  weld  material.   The  volume  of  a  weld  depends  on 
the  size  and  shape  of  the  weld.   For  a  fillet  weld,  the 
volume  is  based  on  a  weld  cross  section  which  is  triangular. 
For  a  butt  weld,  the  volume  is  based  on  the  cross  section 
of  a  double  vee  butt  weld.   F.ach  of  these  welds  is  the 
standard  weld  used  in  highway  girder  design. 

The  objective  function  for  this  problem  is  similar  to 
the  one  found  in  the  CAD- I  program  developed  by  Goble  and 
DeSantis  [19].   This  is  done  so  that  final  comparisons 
between  rol led -sect  ion  girders  and  welded  plate  girders  can 
be  made.   These  comparisons  can  only  be  accurate  if  the 
objective  functions  for  both  types  of  design  arc  similar  in 
form  and  content. 


\ .  7>        Constraints 
The  problem  constraints  are  the  limitations  placed  on 
various  parts  of  the  design  in  order  to  limit  the  total 
design  space.   There  are  two  types  of  constraints  in  an 
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optimization  problem.   One  type  of  constraint,  referred  to 
as  a  side  constraint,  is  imposed  on  the  design  by  functional 
limitations.   The  second  type,  called  a  behavioral  constraint, 
is  imparted  to  the  design  by  a  specific  specification. 

A  functional,  or  side,  constraint  used  in  the  girder 
design  problem  is  the  total  number  of  rolled  sections 
designed.   Since  the  design  of  a  highway  bridge  girder  is 
usually  limited  to  the  36-  and  33-in.  wide-flange  sections, 
the  limited  number  of  design  possibilities  actually  becomes 
a  side  constraint  for  the  problem. 

The  allowable  width  of  the  cover  plates  is  also  con- 
trolled by  a  functional  constraint.   The  plates  are  limited 
to  sizes  of  whole  inches  and  must  be  smaller  than  the  width 
of  the  rol led -section  flange.   The  reason  for  this  is  to 
make  the  welding  process  more  efficient.   With  this  con- 
straint, the  widths  of  the  cover  plates  are  limited  to 
10,  14,  or  15  in. 

The  behavioral  constraints  for  the  highway  girder 
design  are  specified  in  the  American  Association  of  State 
Highway  Officials  bridge  design  code  [20].   The  AASHO  code 
specifies  all  of  the  design  criteria,  which  must  be  converted 
into  the  problem  constraints,  for  the  design  of  a  rolled- 
section  girder.   These  behavioral  constraints  are  divided 
into  the  following  four  groups: 

1  .   Static  stress  constraints, 
2 .   Fatigue  stress  constraints, 
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3.  Cover  plate  size  constraints, 

4.  Secondary  design  constraints. 

The  static  stress  constraints  arc  based  on  the  required 
allowable  stresses  with  no  regard  to  fatigue  considerations. 
There  are  actually  five  different  static  stress  conditions 
which  must  be  satisfied  before  a  design  can  be  considered 
acceptable.   Four  of  the  conditions  consider  the  bending 
stresses  in  the  girder  and  the  fifth  considers  the  shear 
stress  in  the  girder  web.   The  five  constraints  are: 

1.  The  allowable  flexural  tension  stress  in 

the  extreme  fibers  of  rolled  shapes  or  built- 
up  sections  is  0.55  F   where  F   is  the  yield 

v        v        - 

point  stress . 

2.  The  allowable  flexural  compression  stress 
in  the  extreme  fibers  of  rolled  shapes  and 
built-up  sections  is  0.55  F   when  the  com- 

1  y 

pression  flange  is  continuously  supported 
by  being  embedded  in  concrete. 

3.  The  allowable  flexural  compression  stress  in 
the  extreme  fibers  of  rolled  shapes  or  built- 
up  sections  is  determined  by  Equation  4.4: 

3F 


0.55  F 


1.0 


via)21 


TT  2  F.  {  b 


4.4 


This  constraint  reduces  the  allowable  com- 
pressive stress  in  order  to  prevent  the 
buckling  of  the  compressive  flange. 
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4.  The  allowable  compression  stress  in  the 

extreme  fiber-s  of  the  concrete  slab  is 

0.4  f'   where  f '   is  the  28-day  concrete 
c         c 

strength . 

5.  The  allowable  shear  stress  in  the  girder 
web  is  0.35  F  . 

y . 

The  fatigue  stress  constraints  reduce  the  allowable 
stresses  because  of  the  fluctuations  in  stress  caused  by  the 
repeated  applications  of  the  liv.e  load.   The  fatigue  stress 
conditions  are  governed  by  Equations  3.6  and  3.7.   The 
fatigue  theory  used  in  the  design  program  is  described  in 
Section  3.2. 

The  maximum  permissible  thickness  of  a  cover  plate  is 
governed  by  Equation  4.5: 


t       =  1.5  tr  4.5 

cp  f 

*max 

Equation  4.6  governs  the  length  of  the  cover  plate: 


L       =  2  D  +  3  4.6 

cp  • 
1  mm 

where   D  is  the  depth  of  the  rolled  section 
i  n  feet . 
Equations  4.5  and  4.6  set  the  requirements  for  a  cover  plate 
on  a  given  rolled  section  and  make  up  the  cover  plate  size 
constraints  . 

The  final  group  of  design  constraints  control  the 
secondary  details  of  the  girder  design.   The  details  which 
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are  controlled  by  this  group  arc  the  shear  connector  spacings, 
the  bearing  stiffener  design,  and  the  weld  design. 

If  the  design  is  composite,  then  the  following  four 
constraints  control  the  spacing  of  the  shear  connectors: 

1.  The  shear  connectors  are  spaced  so  that  the 
horizontal  shear  is  transferred  from  the  slab 
to  the  steel  girder.   This  spacing  is  designed 
according  to  fatigue  stress  conditions. 

2.  The  design  of  the  shear  connector  spacing  is 
checked  for  ultimate  strength. 

3.  Extra  connectors  must  be  provided  at  the 
points  of  contraf lexure . 

4.  The  maximum  permissible  spacing  is  24  in. 
The  equations  and  the  theory  used  in  the  design  of    these 
shear  connector  constraints  can  be  found  in  Section  5.3. 

The  design  of  the  bearing  stiffeners  is  controlled  by 
the  following  constraints: 

1.  Stiffeners  are  required  at  a  support  if  the 
shear  stress  in  the  web  at  that  support  is 
greater  than  75  percent  of  the  allowable  shear 
for  gi rder  webs . 

2.  The  minimum  thickness  for  a  stiffener  is 
given  by 


m  m 


bi  1  /  ^__ 
12  V  53(100 


I.  ~ 


where  b'  is  the  width  of  the  stiffener. 
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The  bearing  stress  between  the  stiffener 
and  the  rolled-section  flanges  must  be  less 
than  0.8  0  F  . 

y 

The  compressive  stress  in  the  stiffener  must 
be  less  than  the  allowable  column  stress 
determined  from  Equation  4.8: 

'0.75L' 1 2 


0.55F 


y 


1.  25 


1.0 


F 


4tt2E 


4.8 


where   L'/r  is  the  slenderness  ratio  of 

the  stiffener  and  a  section  of  the 
web  which  is  18tw  inches  wide. 
The  final  secondary  design  constraint  is  used  to  con- 
trol the  weld  design.   The  allowable  static  shear  stress  in 
the  weld  material  (Fv)  is  specified  by  Equation  4.9: 


If   F   <  36  ksi,  F   =  12.4  ksi 

y  —        '  v 

If   F   >  56  ksi,  F   =  14.7  ksi 

y  v 


4.9 


The  allowable  shear  stress  for  the  weld  material  must  be 
reduced  for  fatigue  by  using  Equation  3.6  and  the  constants 
under  category  G  in  Table  3.2. 
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CHAPTER  V 

COMPUTER  PROGRAM 

5 . 1    Description  of  the  MA I N  Program 
The  computer  program  for  the  design  of  rol led- section 
girders  consists  of  two  parts.   The  analysis  portion  is 
used  to  determine  the  design  conditions  for  the  problem  and 
the  design  portion  is  used  to  determine  the  optimum  combina- 
tion of  rolled  section  and  cover  plates.   The  entire  program 
is  composed  of  a  MAIN  program  and  thirty- four  subroutines. 

The  purpose  of  the  MAIN  program  is  to  call  on  the 
subroutines  in  the  correct  sequence  to  produce  the  optimum 
design  for  the  particular  girder.   The  following  description 
is  a  step-by-step  summary  of  the  MAIN  program.   A  flowchart 
for  this  program  is  presented  in  Figure  5.1. 

Step  1 
The  rolled  section  information  required  for  the  design 
program  is  read  from  data  cards.   All  the  information  for 
the  eighteen  56-  and  55- in.  wide- flange  sections  is  already 
on  cards  supplied  with  the  program. 

Step  2 
The  READIN  subroutine  is  called  to  read  all  of  the  data 
for  the  design  problem.   A  detailed  description  of  READIN  is 
presented  in  \ppcndix  C. 
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G> 


0- 


Read  the  rolled  section  information 


I 


(2) 


Read  all  the  problem 
data  using  READIN 


f  (5) 


Determine  F 


j  (4) 


Analyze  the  structure 
using  ANAL 


I 


(5) 


CMIN  =  999999 
WMIN  =  999  999 


I 


(6) 


Design  the  uniform  section  grrder, 
I,  with  the  UNIFST  subroutine. 


6 


Figure  5.1:   Flowchart  for  the  MAIN  Program 


"1 


True 


False 


Design  the  rolled  section  J,  using 
the  sequence  of  steps  described  in 
Section  3 . 4 


(8) 

Place  the  optimum  design 
working  arrays 

in  the 

I 


Use  SFMIC  to  determine  the  section 
properties 


t  C£j 


Determine  the  new  girder  weight 


Reanalyze  the  structure 


i;  igure  5.1:   (Cont '  d 
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No 


Print  new  design 
conditions 


t  (11) 


Recalculate  the  stiffeners  and  the 
shear  connector  spacing 


I 


Check  all  the  stresses  using  STRCHK 


Y  e  s 


No 


Print  the  optimum 
design  information 


<D 


igure  5.1:   (Cont ' d . ) 
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Figure    5.1:       (Cont'd. J 
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Determine  the  total  moment  at  the 
critical  section 


Calculate  the  percentage,  PI 


o 


No 


Yes 


Yes 


<5 


i 


ULLL 


Completely  redesign  the  rolled 
section  used  in  the  optimum  design 


I 


(16' 


Calculate  the  deflections  using 
the  DEFLEC  subroutine 


o 


w  (17; 


Print  the  final 
design  information 


figure  5.1:   (Cont'd. ) 
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Step  3 
The  minimum  tensile  strength  used  in  the  fatigue  calcu- 
lations is  determined  from  the  yield  point  stress  of  the 
material  used  in  the  problem.   This  tensile  strength  is  a 
measure  of  the  minimum  ultimate  strength  of  the  steel. 

Step  4 
The  ANAL  subroutine  is  called  to  analyze  the  given 
problem  and  to  determine  all  of  the  design  reactions,  the 
design  moments,  and  the  design  shears.   These  conditions 
are  used  to  determine  the  optimum  girder  for  the  problem. 

Step  5 
The  minimum  weight  and  the  minimum  cost  variables  are 
initialized  at  very  high  values. 

Step  6 
The  UNIFST  subroutine  is  used  to  design  the  uniform 
section  girder.   This  girder  is  the  lightest  rolled  section 
which  satisfies  all  of  the  design  constraints  without  the 
use  of  cover  plates.   This  section  then  becomes  the  initial 
design. 

:  tep  7 
A  design  loop,  which  completely  determines  the  girders 
composed  of  the  sections  lighter  than  the  one  used  as  a 
uniform  section,  is  begun.   This  loop  is  described  in 
Section  3.4  and  a  flowchart  is  presented  in  Figure  3.3. 
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Step  8 
After  the  design  loop  is  completed,  the  optimum  design 
for  the  cycle  is  stored  in  the  optimum  storage  arrays.   This 
design  is  then  placed  into  the  working  arrays  of  the  program 
and  all  of  the  section  properties  are  computed  using  the 
SFMIC  subroutine. 

Step  9 
The  new  unit  dead  weight  of  the  girder  is  determined 
and  the  problem  is  analyzed  again,  using  the  new  distribu- 
tion of  the  moments  of  inertia.   If  the  problem  is  a  con- 
tinuous girder,  then  the  ANAL  subroutine  is  called  to  do  the 
complete  reanalysis.   If  the  problem  is  a  simple- span 
girder,  the  dead-load  design  conditions  are  changed  directly 
in  the  MAIN  program.   This  is  the  only  change  required 
because  the  moments  in  a  simple-span  girder  are  independent 
of  the  cross  -  section  properties. 

Step  10 
If  IPT10  is  greater  than  0,  all  of  the  design  conditions 
are  printed  out . 

Step  11 
Using  the  new,  design  conditions,  the  stiffeners  and 
connector  spacing  for  the  optimum  design  are  recalculated. 
All  of  the  stresses  in  the  girder  arc  determined  bv  the 
STRCliK  subroutine  and  a  new  value  of  the  objective  function 
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is  determined.   All  of  the  information  regarding  this  design 
is  then  printed  if  IPT10  is  greater  than  0. 

Step  12 
If  the  design  cycle  parameter,  IPT7,  equals  2  and  the 
maximum  number  of  cycles  has  not  been  readied,  the  counter 
is  increased  by  one  and  control  is  returned  to  step  5. 
Steps  5  thru  11  are  repeated  until  the  specified  maximum 
number  of  cycles  is  reached.   Control  is  now  passed  to 
step  15. 

Step  13 
If  IPT7  equals  1,  the  total  moment  at  the  critical 
section  is  determined  and  the  percent  change  from  the 
previous  analysis  is  computed.   If  this  percentage  is  less 
than  the  specified  percentage,  control  is  passed  to  step  15. 
Otherwise,  control  is  returned  to  step  5,  where  steps  5  thru 
11  are  repeated  until  the  percentage  is  less  than  that 
specified . 

Step  14 
If  IPT7  equals  0,  no  redesign  cycles  are  required  and 
control  is  passed  to  step  15. 

Step  15 
Using  the  final  conditions,  the  rolled  section  used  in 
the  optimum  design  is  completely  recalculated  using  the 
sequence  of  steps  described  in  Section  3.4. 
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Step  16 

The  required  deflections  are  calculated,  using, the 

moment  of  inertia  distribution  of  the  final  design.   These 

deflections  are  calculated  by  the  DEFLEC  subroutine. 

Step  17 
All  of  the  information  pertaining  to  the  final  design 
is  printed  in  the  output.   This  includes  the  section  proper- 
ties, the  required  cover  plates,  the  required  stiff eners , 
the  shear  connector  spacing  and  all  of  the  deflections.   The 
cost  or  weight  of  the  design  is  also  presented,  depending  on 
the  value  of  1PT6. 

Step  18 
Control  is  returned  to  Step  2  where  the  program  is 
either  terminated  or  the  data  for  a  new  design  problem  are 
read . 


5 . 2    Sample  Problems 

In  order  to  show  the  usefulness  of  the  computer  program, 
two  sample  girders  have  been  designed.   Both  problems  are 
presented  in  the  U.  S.  Steel  Manual  [31].   The  costs  of  the 
manual  design  and  program  design  are  compared  below. 

The  first  problem  is  a  two-span  girder  having  equal 
span  lengths  of  70  ft.   The  design  is  based  on  a  composite 
section  having  the  properties  shown  in  Table  5.1.   The  design 
loading  is  the  AASHO  HS20  loading. 
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Tabic  5.1    Composite  Section  Properties 


Property 

Value 

SLABWD 

84.  00  in 

SLABTI1 

7.00  in 

FPC 

300  0.00  psi 

MR 

8.0  0 

AS 

4.34  sq  in 

HAUNCH 

1.87  5  in 

Table  5.2    Unit  Costs 


Item 


unit 


Cost  in 
Dollars 


Rolled  Section 

Cover  Plates 

Sti  f f eners 

Connectors 

Connector  Installation 

Fixed  Weld  Cost 

Weld  Material  Cost 

Weld  Length  Cost 


lb 

lb 

lb 

lb 

ea 

ea 

cu  in 

i  n 

0.15 
0.12 
0.12 
0.00 
0.75 
20.00 
2.00 
0.0 
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The  optimization  of   the  girder  is  based  on  the  total 
cost.   The  unit  costs  used  in  the  optimization  procedure 
are  shown  in  Table  5.2 

The  final  design  produced  by  the  computer  program  is 
shown  in  Figure  5.2.   The  itemized  costs  and  weights  of 
both  the  computer  design  and  the  manual  design  are  given 
i  n  Table  5.3. 


Table  5.5:   Costs  and  Weights  for  Problem  No.  1 


Item 

U.  S.  Steel  Manual 

Computer  Program 

Weight 

Cost 

Weight 

Cost 

Rolled  Section 
Cover  Plates 
St  if f eners 
Connectors 
Welds 

Totals 

18900.  0 
2361.6 

202.2 
21465.8 

$2855. 00 
$283.40 

$193.  50 
$628.29 

$3940. 19 

22400.0 
243.1 

181.  0 
22824.1 

$3560.00 
$28.16 

$174.00 
$206.90 

$3769.06 

Table  5.5  shows  that  the  optimization  program  produced 
a  girder  which  is  approximately  6.0  percent  heavier  than  the 
girder  found  in  the  manual.   The  optimum  design,  however, 
costs  5.0  percent  less  than  the  manual  design.   Most  of  the 
cost  reduction  can  be  found  in  the  cover  plate  cost.   The 
total  cost  of  the  plates,  including  both  material  and 
fabrication  costs,  for  the  optimum  girder  is  $255.06.   The 
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same  costs  for  the  manual  design  give  a  total  of  $911.69  or 
approximately  four  times  as  much  as  the  program  solution. 
This  cost  difference  offsets  the  increased  cost  of  the 
rolled  section  used  in  the  optimum  design. 

The  second  example  is  a  four-span  composite  girder. 
The  outer  spans  arc  70  ft  and  the  interior  spans  are  90  ft. 
The  composite  section  properties,  the  loading  conditions  and 
the  cost  parameters  are  the  same  as  for  the  first  problem. 
The  program  design  and  the  itemized  costs  and  weights  are 
shown  in  Figure  5.3  and  Table  5.4  respectively. 

Table  5.4:   Costs  and  Weights  for  Problem  No.  2 


Item 

U.  S.  Steel  Manual 

Computer  Program 

Weight 

Cost 

Weight 

Cost 

Rolled  Section 
Cover  Plate 
Stiff eners 
Connectors 
Welds 

Totals 

43200.0 
9855.0 

4  7  5.0 
53520.0 

$6480.00 
$1182.60 

$458.00 
$1598.44 

$9719.04 

51200.0 
3509.6 

548.0 
5  5  0  5  7.6 

$7670.00 
$420.93 

$536.50 
$1251.68 

$9679.11 

Table  5.4  shows  that  the  computer  solution  is  approxi- 
mately 2.9  percent  heavier  and  onlv  0.5  percent  cheaper  than 
the  manual  design.   These  figures  show  that  there  is  a 
slight  saving  in  the  computer  solution.   The  small  amount 
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indicates  that  the  manual  solution  approaches,  for  these 
cost  figures,  the  optimum  design.   The  similar  costs  show 
that  the  computer  not  only  approaches  the  optimum  solution, 
but  also  produces  a  practical  design. 

The  computer  program  is  approximately  4000  cards  in 
length  and  was  developed  using  the  CDC  6500  computer.   The 
compilation  time  for  the  program  is  about  48  seconds.   The 
execution  times  for  the  two  sample  problems  were  60  and  40 
seconds  respectively. 

The  10  x  1/4  plates,  shown  in  Figure  5.3,  violate 
width-to-thickness  ratios  of  the  AASHO  code.   The  program 
is  being  modified  such  that  the  cover  plates  conform  to  the 
code  requirements . 
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CHAPTER  VI 
SUMMARY  AND  CONCLUSIONS 

6  . 1    Summary 

The  cost  optimization  of  a  highway  bridge  girder, 
composed  of  a  rolled  section  acting  compos itely  with  the 
bridge  deck,  has  been  formulated  and  programmed.   An 
objective  function,  based  on  the  material  and  the  fabrica- 
tion costs,  has  been  developed  and  used  in  the  examples. 
The  design  of  the  girder  is  controlled  by  the  1969  AASHO 
code . 

Due  to  the  complex  loadings  required  on  a  highway 
bridge,  the  method  of  influence  lines  is  used  in  the 
analysis.   The  analysis  portion  of  the  program  produces 
the  design  conditions  for  the  problem,  which  are  used  to 
design  and  optimize  the  required  girder. 

Since  the  number  of  possible  designs  is  fairly  limited, 
the  program  uses  the  exhaustive  search  technique.   The 
techniques  of  interval  halving  and  dvnamic  programming  are 
also  used  within  the  design  program.   Interval  halving  is 
used  to  determine  the  thicknesses  of  the  cover  plates  for 
the  subelements  and  dynamic  programming  is  used  to  determine 
the  optimum  location  of  the  cover  plate  splices. 
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The  computer  program  developed  will  solve  a  simple-span 

problem  as  well  as  a  continuous  problem  of  up  to  four  spans. 
The  girder  can  be  either  composite  or  noncompositc  and  is 
designed  for  static  and  fatigue  loads.   All  the  design 
details  are  developed  and  then  listed  in  the  output.   No 
attempt  is  made  to  design  the  girder  for  deflection,  but  the 
total  deflections  for  the  final  design  are  tabulated.   The 
engineer  has  the  option  to  optimize  the  design  based  on 
either  the  total  weight  or  the  total  cost. 

6 . 2    Recommendations  for  Further  Research 

1.  Beams  with  different  steel  strengths  may  be 
cons  idered. 

2.  Varying  the  rolled  section  at  points  of  zero 
dead- load  moment  may  be  considered. 

3.  The  effect  of  field  splices  on  the  optimiza- 
tion problem  could  be  established. 

4.  A  more  accurate  objective  function  could  be 
developed,  provided  that  cost  information 
could  be  obtained  from  steel  fabricators. 

5.  A  complete  optimization  study  of  rolled- 
section  girders  and  plate  girders  could  be 
made  using  the  program  in  conjunction  with 
G  ,D-I  . 

6.  The  effect  of  the  girder  spacing  on  the 
optimum  solution  could  be  determined. 
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The  effect  of  composite  action  in  the  negative 
moment  region  on  the  optimum  solution  could  be 
studied . 

The  effect  of  a  variable  slab  thickness  on  the 
problem  optimization  could  be  determined. 
An  optimization  of  the  bridge  deck  and  stringer 
portion  of  a  highway  bridge  could  be  developed. 
An  optimization  study  of  the  various  types  of 
highway  overpass  structures,  such  as  plate 
girders,  rolled- sec tion  girders,  reinforced 
and  prestressed  concrete  girders  and  concrete 
slab  bridges,  could  be  made  with  the  ultimate 
outlook  on  optimizing  the  entire  bridge 
design  pro j ect . 
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APPENDIX  A 
SUBROUTINE  DESCRIPTIONS 

The  subroutines  used  in  the  program  are  listed  and 
described  in  the  following  list.   The  listing  is  presented 
in  alphabetical  order. 

ALLOW   -  returns  the  allowable  bending  stress 

for  any  location  along  the  top  and  the 

bottom  of  the  girder,  disregarding  all 

fatigue  conditions. 
ANAL    -  determines  the  design  conditions 

(reactions,  moments  and  shears)  for 

a  particular  problem. 
BASFAT  -  returns  the  allowable  base  metal  fatigue 

'stress  at  a  given  point  on  the  girder. 
CONDES  -  designs  the  shear  connectors  for  a  given 

girder . 
CPDES   -  designs  the  cover  plates  for  a  given 

rolled  section. 
CPWDES  -  designs  the  required  end  and  seal  welds 

for  each  cover  plate  required  on  the 

girder . 
DEFLEC  -  determines  the  required  deflections  for 

the  final  girder  design. 
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EQSET   -  calculates  the  simple  beam  deflections, 
at  the  interior  support  points,  which 
are  used  in  calculating  the  influence 
lines  for  the  reactions. 

FLOPT   -  uses  the  optimization  method  of  dynamic 
programming  to  determine  the  final  cover 
plate  arrangement  along  a  girder. 

ICALC   -  determines  the  moments  of  inertia  and 

the  distances  to  the  lower  extreme  fibers 
for  any  rolled  section  and  cover  plate 
combination.   It  is  also  used  to  find  the 
similar  properties  for  the  composite 
sections  by  treating  the  slab  as  a  top 
cover  plate. 

ILINT   -  returns  the  interpolated  value  of  any 
point  along  an  influence  line  using 
straight  line  interpolation. 

ILPROP  -  determines  the  various  required  proper- 
ties of  a  given  influence  line  including 
the  crossover  points,  the  positive  and 
negative  areas  and  the  locations  of  the 
maximum  and  minimum  ordinates. 

IMPACT  -  determines  the  impact  factor  for  any 
particular  analysis  function. 

KEEP    -  decides  whether  or  not  to  retain  the 
given  design  as  the  present  optimum. 
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LOAD    -  places  all  the  required  loads  on  the 

influence  line  in  such  a  manner  as  to 

return  the  maximum  loading  conditions 

required. 
MAIN    -  calls  on  the  remaining  subroutines  in 

the  correct  order  to  produce  the  optimum 

design. 
OBJECT  -  contains  the  function  being  minimized 

during  the  optimization  of  the  design 

problem . 
PRCON   -  prints  the  connector  design  for  a  given 

girder . 
PRCOV   -  prints  the  cover  plate  information  for  a 

girder . 
PRDEF   -  prints  the  deflections  for  the  final 

girder  design . 
PRSEC   -  prints  the  section  table  for  a  given 

girder . 
PRSTIF  -  prints  the  bearing  stiffener  information 

for  a  given  girder  design. 
PRSTR   -  prints  the  stress  table  for  a  given 

g  i  rder  des  i  gn . 
REAC    -  determines  the  reaction  influence  lines 

for  a  continuous  girder  using  the  moment- 
area  theorems . 


READIN  -  reads  in  all  the  pertinent  data  required 
for  a  particular  design  problem. 

SFMIC   -  calculates  the  properties  for  each  sub- 
element  along  the  girder. 

SORT    -  places  the  required  cover  plate  thick- 
nesses in  ascending  order  so  they  can 
be  used  in  the  FLOPT  subroutine. 

STIFF   -  completely  designs  the  bearing  stiffeners 
for  a  given  girder  if  they  are  required. 

STRCHK  -  determines  the  stresses  on  a  given  girder 
and  checks  them  against  the  allowable 
stresses . 

STRMAX  -  determines  the  maximum  design  stress 
conditions . 

TERMIN  -  determines  the  locations  of  the  ends  of 
the  required  cover  plates  based  on  the 
fatigue  stress  conditions. 

'HIDES   -  calculates  the  required  cover  plate 

thicknesses  for  each  subelement  along 
the  g i  rder . 

TOPT    -  determines  the  required  minimum  weld 
thickness  Lor  a  fillet  w  e 1 d . 

UNIFST  -  designs  the  girder  which  has  a  uniform 

steel  section  composed  of  a  rolled,  ki de- 
fiance section  without  cover  plates. 


APPENDIX    B 


95 


APPENDIX  B 
SELECTED  PROGRAM  NOMENCLATURE 

The  following  alphabetic  list  contains  some  of  the 
more  pertinent  variables  used  in  the  program.   The  list 
basically  contains  all  the  variables  which  are  interchanged 
between  the  program  subroutines.   The  numbers  in  the  paren- 
theses immediately  following  the  \rariable  name  give  the 
maximum  dimensions  of  that  array  within  the  program. 

ANEG  -  the  negative  area  under  a  given 

influence  line. 
APOS  -  the  positive  area  under  a  given 

influence  line. 
BEAR  (5, 3  J     -  the  bearing  stiffener  array. 

(1.1)  -  width  of  stiffener  at  support  I. 

(1.2)  -  thickness  of  stiffener  at  support  I 

(1.3)  -  size  of  connecting  weld  for 

stiffener  at  su^ort  I. 

CC    (5,5)     -  the  coefficient  array  for  the  reac- 
tion influence  line  calculations. 

CH    (80)      -  the  distance  from  the  neutral  axis 

to  the  lower  extreme  fiber  within 
a  sub  el  orient  for  the  high  modulus 
composite  section. 
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CL    (80)      -  the  distance  from  the  neutral  axis 

to  the  lower  extreme  fiber  within 

a  subelement  for  the  low  modulus 

composite  section. 
CMIN  -  the  cost  of  the  present  optimum 

design. 
CONSP  (1 6 , 3 , 4 )  -  the  shear  connector  design  array. 

(J,  1,1)   -  spacing  of  the  J  group  of  connectors 

in  span  1 . 
(J,2,I)   -  the  starting  coordinate  of  the  J 

group  of  connectors  in  span  I . 
(J, 3, I)   -  the  final  coordinate  of  the  J  group 

of  connectors  in  span  I . 
COOR  (81)      -  the  coordinates  of  the  analysis 

points . 
COST  (9)       -  the  unit  cost  array. 

(1)  -  cost  of    the  rolled  sections. 

(2)  -  cost  of  the  cover  plates. 

(3)  -  cost  of  the  bearing  stiffeners. 

(4)  -  cost  of  the  shear  connectors. 

(5)  -  cost  of  the  connector  installation. 

( 6 )  -  fixed  cost  of  a  weld. 

( 7  J       -  cost  of  the  weld  material. 
(8)         cost  of  the  weld  Length. 
COVPL (12 , 6 , 2)  -  the  cover  plate  design  array. 

(J,  1,1)   -  the  starting  coordinate  of  Jth  plate 
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(J  ,2,1) 
(J  ,  3 ,  I  ) 
(J  ,4,D 
(J,  5,  I) 


CS    (80) 


DAF   (18) 


the  final  coordinate  of  Jth  plate. 

the  thickness  of  Jth  plate. 

the  width  of  the  Jth  plate. 

the  size  of  the  end  welds  for  the 

Jth  plate. 
(J  ,  6 , I )   -  the  size  of  the  seal  weld  for  the 

Jth  plate. 
Note:   For  the  COVPL  array,  I  equals  1  for  the 
cover  plates  on  the  top  of  the  girder 
and  I  equals  2  for  the  cover  plates  on 
the  bottom  of  the  girder. 

the  distance  from  the  neutral  axis 

to  the  lower  extreme  fiber  within 

a  subelement  for  the  steel  section 

alone . 

the  depth  divided  by  the  flange 

area  for  the  rolled  sections  used 

in  the  design  program. 

the  depth  of  the  rolled  sections 

u  s  e  d  i  n  t  h  e  d  e  s  i  g  n  p  r  o  g  r  am  . 

the  design  girder  deflections. 

the  deflection  at  I  OFF  ( I )  due  to 

the  dead  load . 

the  deflection  at  mEF(I)  due  to 

the  positive  live  load. 


DEPTH  (1  S) 

DESDEF (2  0,4) 
(1,1) 

(1,2) 
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(1.3)  -  the  deflection  at  IDEF(I)  due  to 

the  superimposed  dead  load. 

(1.4)  the  deflection  at  IDEF(I)  due  to 
the  negative  live  load. 

DESMOM (81 , 4 )   -  the  design  moment  array. 

(1.1)  -  the  moment  at  I  due  to  the  dead 

load  . 

(1.2)  -  the  moment  at  f  due  to  the  positive 

live  load . 

(1.3)  -  the  moment  at  I  due  to  the  super- 

imposed dead  load. 

(1.4)  -  the  moment  at  I  due  to  the  negative 

1 ive  load . 
DESREA(5,4)    -  the  design  reaction  array. 

(1.1)  -  the  reaction  at  support  t  due  to 

the  dead  load. 

(1.2)  -  the  reaction  at  support  I  due  to 

the  positive  live  load. 

(1.3)  -  the  reaction  at  support  1  due  to 

the  superimposed  dead  load. 

(1.4)  -  the  reaction  at  support  I  due  to 

the  negative  live  load. 
DESSH  (162,4)  -  the  design  shear  array. 

(1,1)    -  the  shear  at  point  I  due  to  the 
d  c  a  d  1  o  a  d  . 
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(1,2 J         -    the   shear   at   point    I    clue    to    the 

positive  live  load. 
(1,3)    -  the  shear  at  point  I  due  to  the 

superimposed  dead  load. 
(I,4J    -  the  shear  at  point  I  due  to  the 

negative  live  load. 
Note:   For  analysis  point  I,  J.  equals  21-1 
just  to  the  right  of  I  and  J  equals 
21  -  2  just  to  the  left  of  I. 
DIA  -  the  diameter  of  the  shear  connector 

used  in  the  program. 
DING  -  the  impact  length  for  a  given 

influence  line . 
DIST  -  the  distance  a  cover  plate  is  cut 

off  based  on  fatigue  stress  con- 
siderations . 
EM  -  the  modulus  of  elasticity  which  is 

set  at  29,000  ksi  in  the  program. 
EPS    (4)      -  the  number  of  elements  per  span. 
FC  -  the  allowable  concrete  compressive 

stress  . 
FLTH   (IS)     -  the  thickness  of  the  flange  for 

the  rolled  sections  used  in  the 


i  i  w in   (is 


d  e  s  l  g  n   p  r  o  g  r  a  n  . 

the   width    of    the    flange    for    the 

rolled    sections    used    in    the    design 


pre  rar 


LOO 

FPC  -  the  28-day  strength  of  the  concrete 

I-'U  -  the  minimum  tensile  strength  of 

the  steel . 

FV  -  the  allowable  weld  metal  stress. 

GDL  -  the  weight  of 'the  steel  girder. 

11AUNCH        -  the  value  of  the  concrete  haunch 

for  the  composite  design. 

1IT  -  the  calculated  impact  factor. 

IDEF   (30)     -  the  analysis  points  at  which  the 

deflections  are  calculated. 

IHCON  (80)     -  the  moment  of  inertia  of  a  sub- 
element  for  the  high  modulus 
composite  section. 

ILCON  (80)     -  the  moment  of  inertia  of  a  sub- 
element  for  the  low  modulus 
composite  section. 

IOSP   (5)      -  the  index  of  the  support  points. 

IPT1  -  the  subelement  input  option. 

IPT2  -  the  loading  type  input  option. 

IPT3  -  the  initial  cross  section  input 

opt  ion . 

IPT4  -  the  units  input  option. 

IPT5  -  the  loading  components  or  designa- 

tion input  option. 

IPT6  -  the  optimization  input  option. 

[PT7  -  the  design  cycle  input  option. 
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I  P'1'8  the  materia]  input  option. 

[PT9  -  the  deflection  input  option. 

1IT10  -  the  output  option. 

[SAVE  (80)     -  the  rolled  section  number  for  each 

sub clement  of  the  present  optimum 
des  ign  g i  rder . 

1ST  -  the  identification  number  of  the 

rolled  section  being  designed. 

ISTEEL(80)     -  the  moment  of  inertia  of  a  sub- 
element  for  the  steel  section. 

IX     (18)       the  moment  of  inertia  about  the 

x-x  axis  for  the  rolled  sections 
used  in  the  design  program. 

IV     (18)     -  the  moment  of  inertia  about  the 

y-y  axis  for  the  rolled  sections 
used  in  the  program. 

JPT3  -  the  composite  section  variable. 

Note:   If  the  design  is  composite  then  JPT3 
is  less  than  9  and  if  the  design  is 
noncomposite  then  JPT3  is  greater 
than  9. 

KOP    fin)     -  the  location  oi'    the  crossover 

points  for  a  given  influence  line. 

Li)  IS  -  the  live  load  designation. 

LENGTH  the  total  length  of  the  girder. 
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LPCM   (10) 


MAXOL 


M  I  NOL 


MN 


MR 

NA 

NAME   (18) 

NCS 

NCY 
NCYC 

NDEF 

NE 

NOS    (8  0) 

N.S 

NSECT  (18) 


the  boundary  points  of  the  negative 

moment  regions  along  the  girder. 

the  point  which  is  the  maximum 

ordinate  of  a  given  influence  line. 

the  point  which  is  the  minimum 

ordinate  of  a  given  influence  line. 

the  number  of  spans  along  the 

design  problem  girder. 

the  modular  ratio. 

the  total  number  of  analysis 

points  along  the  girder. 

the  AiSC  name  of  the  rolled  sections 

used  in  the  design  program. 

the  number  of  shear  connectors  per 

row. 

the  design  cycle  counter. 

the  maximum  number  of  design 

cycles  required. 

the  total  number  of  points  at 

which  the  deflections  are  required. 

the  total  number  of  subclements . 

the  NSECT  number  of  the  rolled 

section  in  each  subelement. 

the  total  number  of  supports. 

the  identification  number  of  the 

rolled  sections  used  in  the  design 

program . 


lo: 


01 

0  I  L 
PERCEN 

PLATH  (80,2) 

(I,D-  • 

(1,2) 

PLL 
PLTSAV(80,2) 

PLWID  (80,2) 
(1,1) 

(1,2) 

PL WS AY (8 0,2) 

PRO]1   (  5  ) 

I  1  ) 
(2) 
(3  l 
RIL    I 5  ,S]  ,3) 


the  Impact  option. 

the  interstate  loading  option. 

the  allowable  percentage  change 

in  the  critical  design  moment. 

the. plate  thicknesses  for  each 

subelement . 

the  thickness  of  the  top  cover 

plate  in  subelement  I . 

the  thickness  of  the  bottom  plate 

in  subelement  I . 

the  live  load  factor. 

the  PLATH  array  for  the  current 

optimum  design. 

the  plate  width  for  each  subelement 

the  width  of  the  top  cover  plate 

in  subelement  I . 

the  width  of  the  bottom  cover 

plate  in  subelement  1. 

the  PLWID  array  for  the  current 

optimum  design. 

t h e  steel  properties  u s e d  in  the 

d  e  s  i  g  n  . 

the  allowable  bending  stress. 

the  allowable  shear  stress. 

the  steel  yield  point  stress. 

the  reaction  influence  line  array. 
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(I,J,1)  -  the  reaction  at  support  1,  with 

the  unit  load  at  point  .1,  for  the 

steel  section. 
(I, J, 2)  -  the  reaction  at  support  T,  with 

the  unit  load  at  point  J,  for  the 

low  modulus  composite  section. 
(I, J, 5)  -  the  reaction  at  support  I,  with 

the  unit  load  at  point  .T ,  for  the 

high  modulus  composite  section. 
RATIO         -  the  ratio  of  the  minimum  stress 

to  the  maximum  stress  used  in  the 

fatigue  calculations. 
SAREA  (18)     -  the  area  of  the  rolled  sections 

used  in  the  design  program. 
SAVE   (12,6,2)-  the  cover  plate  array  (COVPL)  for 

the  current  optimum  design. 
SDL  -  the  superimposed  dead  load. 

SLABA         -  the  area  of  the  concrete  slab 

which  acts  compositely  with  the 

steel  section. 
SLABTH        -  the  thickness  o\~    the  concrete  slab. 
SLABWD        -  the  effective  width  of  the  concrete 

slab. 
SLABWT         -  the  weight  of  the  concrete  deck 

which  is  carried  by  the  g i r d e r . 
SLL  -  the  sidewalk  live  load. 
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SMAX  -  the  maximum  stress  at  a  point  used 

to  determine  RATIO. 
SMIN  -  the  minimum  stress  at  a  point  used 

to  determine  RATIO. 
STORE  (10)     -  the  loadings  for  a  given  influence 

line. 
STRESS  (81  ,4,2)-  the  stress  array. 

(1,1,1)  -  the  bottom  steel  stress  at  point  I. 
(1,2,1)  -  the  top  steel  stress  at  point  I. 
(1,3,1)  -  the  concrete  stress  at  point  I. 

(1.4.1)  -  the  shear  stress  at  point  I. 

(1.1.2)  -  the  bottom  stress  indicator. 
(1,2,2)  -  the  top  stress  indicator. 
(1,5,2)  -  the  concrete  stress  indicator. 
(1,4,2)  -  the  shear  stress  indicator. 
Note:   If  the  indicator  is  1.0,  the  stress 

is  greater  than  the  allowable  and 

if  the  indicator  is  0.0,  the  stress 

is  less  than  the  allowable. 

SUBLEN(80)     -  the  subelement  length  array. 

TCOST         -  the  total  cost  of  a  girder  design. 

TDS  -  the  total  design  stress  at  a  point. 

TITLE  (36)     -  the  problem  title  array. 

TRLO   (10)     -  the  truck  load  components. 

(1)  -  the  first  axle  load. 

( 2 )  -  the  second  axle  load. 
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(5)      -  the  third  axle  load. 

(4)  -  the  spacing  between  the  first  and 

second  axles . 

(5)  -  the  minimum  spacing  between  the 

second  and  the  third  axles. 

(6)  -  the  maximum  spacing  between  the 

second  and  third  axles. 

(7)  -  the  superimposed  dead  load. 

(8)  -  the  weight  of  the  steel  girder. 

(9)  -  the  impact  option. 

(10)  -  the  weight  of  the  slab  carried  by 

the  steel  girder. 
TSWID1        -  the  transformed  slab  width  for  the 

low  modulus  composite  section. 
TSWID2         -  the  transformed  slab  width  for  the 

high  modulus  composite  section. 
UNLO   (9)      -  the  uniform  load  components. 

(1)  -  the  equivalent  uniform  live  load. 

(2)  -  the  concentrated  load  for  moment. 

(3)  -  the  concentrated  load  for  shear. 

(4)  -  the  superimposed  dead  load. 

(5)  -  the  weight  of  the  steel  girder. 

( 6 )  -  the  interstate  concentrated  load. 

(7)  -  the  sidewalk  live  load. 
(X)      -  the  impact  option. 

(9)      -  the  weight  of  the  slab  carried  by 
the  steel  gird c  r . 
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IVEBTH  (18)     -  the  thickness  of  the  webs  for  the 

rolled  sections  used  in  the  design 

program . 
IVMIN  -  the  weight  of  the  present  optimum 

design . 
WTOT  -  the  total  weight  of  a  given  girder 

design . 


APPENDIX    C 
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APPENDIX  C 
PROGRAM  DATA 

READ  IN  Subroutine 
The  READIN  subroutine  is  used  to  provide  the  program 
with  all  of  the  necessary  data  to  solve  a  given  problem. 
The  routine  docs  not,  however,  read  the  data  cards  which 
pertain  to  the  section  properties  of  the  wide-flange  rolled 
sections.   These  cards  are  read  by  the  main  program  so  that, 
if  multiple  problems  are  being  solved  with  one  run  of  the 
program,  the  section  properties  are  read  only  once.   The 
order  of  operations  used  in  the  READIN  subroutine  are 
described  below.   After  the  description,  the  method  of  data 
input  is  explained  and  a  few  samples  of  input  data  are 
presented . 

Step  1 
All  of  the  arrays  involved  in  the  subroutine  arc 
zeroed  and  prepared  for  the  given  problem. 

Step  2 
The  title  cards  are  read  and  placed  in  the  TITLE  array. 
Two  cards  must  be  present  for  the  title  description.   If  a 
title  is  not  desired  or  if  the  title  is  only  on  one  card, 
then  blank  cards  must  appear  in  this  position. 
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Step  3 
The  input  parameter  card,  containing  the  input  param- 
eters IPT1  to  IPT10,  is  read.   These  parameters  are  defined 
in  Appendix  B. 

Step  4 
The  subroutine  prints  out  some  specified  heading 
information  if  the  output  parameter,  IPT10,  dictates  this 
operation . 

Step  5 
The  lengths  of  the  individual  spans  are  read  and  used 
to  determine  the  total  length  of  the  girder  and  the  number 
of  supports . 

Step  6 
Depending  on  the  value  of  the  subelement  input  option, 
IPT1,  the  subelement  information  is  either  read  from  cards 
or  determined  by  the  program.   The  quantities  included  in 
the  subelement  information  are  the  number  of  elements  per 
span,  the  analysis  point  coordinates,  the  support  point 
indices,  and  the  lengths  of  the  subelcments.   If  the  design 
is  continuous  and  composite,  the  limits  of  the  negative 
moment  regions  arc  tentatively  defined  as  the  two  elements 
on  either  side  of  the  interior  supports.   These  limits  are 
then  stored  in  the  LPCM  array. 
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Step  7 
The  data  determined  in  steps  5  and  6  are  printed  out 
and  converted  to  inches  if  required.   This  conversion  is 
necessary  because  the  program  operates  in  the  units  of  kips 
and  inches . 

Step  8 
The  steel  table  information  for  the  wide-flange  rolled 
sections  used  in  the  program  is  printed  if  the  output  param- 
eter, IPT10,  indicates  that  this  'information  is  desired. 

Step  9 
If  the  design  is  composite,  then  the  composite  section 
data  are  read  and  printed.   The  composite  action  parameter, 
JPT3 ,  is  calculated.   If  the  design  is  noncomposi te ,  step  9 
is  omitted. 

Step  10 
Depending  on  the  value  of  IPT3,  the  initial  cross 
section  for  each  subelemcnt  along  the  girder  is  either  read 
from  cards  or  determined  by  the  subroutine.   With  the 
sections  specified  for  each  subelemcnt,  SFMIC  is  called  to 
evaluate  the  properties  for  each  clement.   These  properties 
are  then  printed  in  tabular  form  if  required  by  the  output 
parameter . 


Step  11 
The  required  loading  conditions  desired  in  the  problem 
are  read  from  cards.   The  form  of  the  input  depends  on  the 
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value  of  the  loading  parameter,  IPT2,  and  the  value  of  the 
loading  components  parameter,  JPT5.   The  TRLO  and  the  II.M 
arrays  are  converted  to  the  units  of  inches,  if  necessary, 
and  are  then  printed  in  the  output. 

Step  12 
The  properties  of  the  steel  used  in  the  problem  are 
either  read  or  determined  by  the  routine.   If  the  steel 
used  has  the  ASTM  designation  of  A36,  then  the  PROP  array 
is  completely  prepared  by  the  subroutine.   For  all  other 
steels,  the  yield  point  stress  is  read  from  a  data  card  and 
the  properties  are  calculated.   The  value  of  the  modulus  of 
elasticity  for  the  steel  is  set  at  29,000  ksi. 

Step  13 
The  allowable  compression  stress  for  the  concrete  is 
computed  using  the  28-day  concrete  strength  read  in  with 
the  composite  section  data.   All  of  the  steel  and  concrete 
allowable  stress  properties  are  then  listed  in  the  output. 

Step  14 
If  the  optimization  is  based 'on  cost,  the  unit  cost 
data  card  is  read.   This  card  is  not  read  if  the  optimization 
is  based  on  weight.   The  COST  array  is  then  added  to  the 
output . 


Step  15 
The  method  of  terminating  the  design  program  is  either 
read  or  determined  by  the  routine.   The  method  which  is 
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employed  depends  on  the  value  of  the  design  cycle  parameter, 
[PT7.   This  termination  method  is  then  described  in  the 
output. 

Step  16 

The  points  at  which  the  deflections  are  required  are 
either  read  or  determined  by  the  subroutine.   The  form  of 
this  input  depends  on  the  value  of  the  deflection  parameter, 
IPT9.   The  points  at  which  the  deflections  will  be  calculated 
are  then  listed  in  the  output. 

This  completes  the  READ  IN  subroutine.   Once  the  data 
have  been  completely  determined,  the  control  of  the  program 
is  returned  to  the  main  program  and  the  design  and  optimiza- 
tion of  the  desired  problem  is  started. 


L13 


Method  of  Data  Input 
The  method  of  data  input  for  this  program  resembles 
the  method  used  in  the  plate  girder  design  program,  GAD-I. 
It  should  be  noted,  however,  that  although  the  methods  are 
similar,  the  data  cards  cannot  be  interchanged.   The  reason 
is  that  some  of  the  data  have  a  different  meaning  and 
purpose  in  this  program.   The  designer  should  note  these 
differences  and  take  the  appropriate  precautions  when  work- 
ing with  both  programs.   The  following  is  a  description  of 
the  method  of  data  input. 

I.   Card  Set  1  -  Title  Cards 

A.   The  card  set  consists  of  two  cards  with 
the  first  72  columns  used  on  each  card. 
Format  (18A4/18A4) 
II.   Card  Set  2  -  Input  Parameter  Card 

A.   The  card  set  consists  of  one  card  with 
the  following  integer  information: 
1 .   IPT1  -  Subelement  Information 

0  -  10  elements  per  span. 

1  -  Head  in  the  coordinates. 

2  -  Read  in  the  number  of  elements 


per  span. 

1 1'T 2  -  Type  of  Loading  Required 
0  -  Lane,  Interstate,  Sidewalk  and 
Truck  L oad ins . 
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1  -  Lane,  Interstate  and  Sidewalk 

Loading . 

2  -  Truck  Loading. 

IPT3  -  Initial  Girder  Section 

0  -  W53xll8  used  with  composite 

des  ign . 

1  -  Read  the  uniform  section  with 

composite  design. 

2  -  Read  the  section  for  each  sub- 

element  with  composite  design. 

10  -  W35*118  used  with  noncomposite 

des  ign . 

11  -  Read  the  uniform  section  with 

noncomposite  design. 

12  -  Read  the  section  for  each  sub- 

element  with  noncomposite  design 
IPT4  -  Units  of  Data 

0  -  Units  of  inches. 

1  -  Units  of  feet. 

IPT5  -  Loading  Components/Designation 

0  -  Read  in  the  load  factor  and 

the  truck  designation. 

1  -  Read  in  all  the  load  components. 
IPT6  -  Type  of  Optimization 

0  -  Optimize  based  on  cost. 

1  -  Optimize  based  on  weight. 
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[PT7  -  Design  Cycles 

0  -  Use  only  one  cycle  of  design. 

1  -  Base  the  design  on  the  percentage 

change  in  the  critical  moment. 

2  -  Read  in  the  maximum  number  of 

des  i  gn  cycles . 
IPT8  -  Steel  Used  in  the  Design 

0  -  Use  A56  steel . 

1  -  Read  the  required  yield  point 

stress . 
1PT9  -  Deflections  Required 

0  -  Determine  the  deflections  at 

the  span  centerlines. 

1  -  Determine  the  deflections  at 

the  0.2,  0.5,  and  0.8  points  in 
each  span. 

2  -  Read  in  the  analysis  points  at 

which  the  deflections  are 
requi red . 
1PT10  -  Output  Parameter 

0  -  Data  and  the  final  solution. 

1  -  Level  0  plus  each  cycle  solution. 

2  -  Level  1  plus  each  girder  design. 
5  -  Level  2  plus  additional  output. 

4  -  Debugging  output. 
ormat  (1015) 
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III.   Card  Set  3  -  Span  Length  Card 

A.   The  card  set  consists  of  one  card  which 
contains  the  lengths  of  the  spans. 
Format  (4F10.4) 
IV.   Card  Set  4  -  Subelement  Information 

A.  If  IPT1  =  0,  card  set  4  is  not  required. 

B.  If  IPT1  =  1,  the  card  set  consists  of  the 
following  cards: 

1 .  One  card  containing  the  total  number 
of  analysis  points. 

Format  (110) 

2.  The  necessary  number  of  cards,  each 
containing  seven  coordinates,  to  com- 
pletely define  the  total  number  of 
analysis  points. 

Format  (7F10.2) 

C.  If  IPT1  =  2,  read  the  number  of  elements 
for  each  span. 

Format  (4110) 
V.   Card  Set  5  -  Composite  Section  Requirements 

A.  If  IPT3  =  10,  11,  or  12,  card  set  5  is 
not  required. 

B.  If  IPT3  =  0,  1,  or  3,  read  the  following 
composite  section  requirements  from  one 
card  : 

1.   Effective  width  o^    the  slab,  in. 
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2.   Thickness  of  si  ah,  in. 
5.   Modular  Ratio  as  a  floating  point 
number . 

4.  28 -day  concrete  strength,  psi. 

5.  Slab  steel  area,  sq  .  in. 

6.  Value  of  the  haunch,  in. 
Format  (5F1 0  .  2  ,  Fl 0  .  4  J 

VI.   Card  Set  6  -  Initial  Cross  Sections 

A.  If  IPT3  =  0  or  10,  card  set  6  is  not 
requi  red . 

B.  If  IPT5  =  1  or  11,  the  card  set  consists 

of  one  card  containing  the  following  values 

1.  ID  number  of  the  rolled  section. 

2.  Width  of  top  cover  plate,  in. 

3.  Thickness  of  top  cover  plate,  in. 

4.  Width  of  the  bottom  cover  plate,  in. 

5.  Thickness  of  the  bottom  cover  plate,  in 
Format  (I10,4F10.4) 

C.  If  IPT5  =  2  or  12,  the  card  sets  consist 
of  one  card  for  each  subelement.   The 
following  values  are  contained  on  each 
card  : 

1.  IP  number  of  the  rolled  section. 

2 .  W i d t h  of  to p  cover  plate,  in. 

3.  Thickness  of  top  cover  plate,  in. 

4.  Width  oi~    the  bottom  plate,  in. 


1  IS 


5.   Thickness  of    the  bottom  plate,  in. 
Format  ( 1 10  ,4F 10 . 4 ) 
VII.   Card  Set  7  -  Load  Information 

A.   If  IPT5  =  0,  the  card  set  consists  of  one 
card  containing  the  following  values: 

1.  Load  designation- 

a.  Use  one  of  the  five  (5)  standard 
AASHO -designations ,  i.e.,  1110, 
HIS,  JI20,  HS15  or  IIS20. 

b.  This  input  must  be  right 
justified  in  the  first  4  columns. 

2.  Load  factor  in  decimal  form. 

3.  Weight  of  the  slab  carried  by  the 
girder . 

4.  The  superimposed  dead  load  acting  on 
the  low  modulus  concrete  section. 

5.  Dead  weight  of  the  girder. 

6.  Sidewalk  live  load. 
7  .   Impact  option . 

a.  Use  0.0  if  impact  is  included. 

b.  Use  1.0  if  impact  is  to  be 
excluded  from  the  design. 

S.   Interstate  , loading  option. 

a.   Use  0.0  if  loading  is  to  he 
cons  i do  red . 


1  1  9 


b.   Use  1.0  if  loading  is  excluded. 

Format  (A4  ,  F6  .  3 ,  6F1  0  .  4  J 
If  [PT5  =  1  and  IPT2  =  0,  the  card  set 
consists  of  two  cards.   The  first  card 
contains  the  truck  loading  components 
(TRLO)  and  the  second  card  contains  the 
uniform  load  components  (UNLO) . 

1 .  Truck  load  components . 

a.  First  axle  load  in  kips. 

b.  Second  axle  load  in  kips. 

c.  Third  axle  load  in  kips. 

d.  Spacing  of  the  first  two  axles. 

e.  Minimum  spacing  between  the 
rear  axles. 

f.  Maximum  spacing  between  the 
rear  axles. 

g.  Superimposed  dead  load. 
h.   Girder  dead  weight. 

i .   Impact  opt  ion . 
i .   Weight  of  slab . 
Format  (10F7.3) 

2.  Uniform  load  components. 

a.  Uniform  lane  load. 

b.  Concentrated  load  for  moment. 

c.  Concentrated  load  for  shear. 

d.  Superimposed  dead  load. 
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c.   Girder  dead  weight. 

f.  Interstate  concentrated  load. 

g.  Sidewalk  live  load, 
h .   Impact  opt  i  on . 

i.   Weight  of  slab. 
Format  (9F8.4) 

C.  If  IPTS  =  1  and  IPT2  =  1,  the  card  set 
consists  of  one  card  containing  the  lane 
loading  components  shown  above. 

D.  If  IPT5  =  1  and  IPT2  =  2,  the  card  set 
consists  of  one  card  containing  the  truck 
loading  components  shown  above. 

VIII.   Card  Set  S  -  Material  Properties 

A.  If  IPTS  =  0,  card  set  8  is  not  required. 

B.  If  IPTS  =  1,  the  card  set  consists  of  one 
card  containing  the  yield  point  stress  of 
the  steel  used  in  the  problem.   Up  to 

4  percent  over stress  can  be  programmed  by 
increasing  the  yield  point  stress  by  the 
desired  overstress  percentage. 
IX.   Card  Set  9  -  Optimization  Information 

A.   If  IPT6  =  0,  the  card  set  consists  of  one 
card  containing  the  following  values: 
1.   Rolled  section  cost,  S/lb. 
2  .      Cover  plate  cost,  $ / 1 b . 
5.   lie;)  ring  stiffener  cost,  $/lb. 


1  2] 


4.  Shear  connector  cost,  $/lb. 

5.  Connector  installation  cost,  $/ea. 

6.  Fixed  weld  cost,  $/weld. 

7.  Weld  material  cost,  $/cu.  in. 

8.  Yield    length  cost,  S/in. 
Format  (9F8.2) 

B.   If  IPT6  =  1,  card  set  9  is  not  required. 
X.   Card  Set  10  -  Design  Cycle  Information 

A.  If  IPT7  =  0,  card  set  10  is  not  required. 

B.  If  IPT7  =  1,  the  card  set  consists  of  one 
card  containing  the  following  values: 

1.  Maximum  number  of  cycles. 

2.  Allowable  percentage  change  in  the 
critical  moment. 

Format  (I10,F10.2) 

C.  If  IPT7  =  2,  the  card  set  consists  of  one 
card  containing  the  number  of  design 
cycles  desired. 

Format  (HO)   - 
XI.   Card  Set  11  -  Deflection  Information 

A.   If  [PT9  =  0,  card  set  11  is  not  required. 
B..  If  IPT9  =  1,  card  set  11  is  not  required. 
C.   If  1PT9  =  2,  the  card  set  consists  of  one 
card  containing  the  following: 
1.   Number  of  points  at  which  the  deflec- 
tions are  required.   (Max.  =  20) 
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2.   Analysis  point  numbers  of  the 
deflection  points. 
Format  (.110,20131 

Note:    When  preparing  the  data,  all  integer 

values  should  be  right  justified  and  all 
floating  point  numbers  should  be  typed 
according  to  their  respective  formats. 
All  floating  point  fields  are  large 
enough  to  include  a  punched  decimal 
point  and  this  is  recommended  to  avoid 
confus  ion . 
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Input  Data  Samples 


Problem  No .  1 

Design  a  girder  with  three  spans,  each  having  a  length 
of  70  ft.   Design  the  girder  for  HS20  loading  with  a  load 
factor  of  one  axle  per  girder.   Also  design  for  impact 
loading  and  interstate  loading.   Use  A36  steel  in  the  design, 
and  base  the  design  on  ten  elements  per  span. 

Design  the  girder  compositely  with  the  deck  slab.   Use 
concrete  with  a  28-day  strength  of  5,000  psi  and  a  modular 
ratio  of  8  for  the  bridge  deck.   The  composite  section 
properties  include  an  effective  slab  width  of  84  in.,  a 
slab  thickness  of  7  in.,  a  slab  steel  area  of  1.94  sq.  in., 
and  a  haunch  depth  of  1.5  in. 

The  optimization  is  based  on  cost.   The  unit  cost  per 
pound  of  steel  for  all  steel  elements  is  10  cents.   The 
fixed  weld  cost  is  $30.00.   The  weld  material  costs  $5.00 
per  cubic  inch  of  material  and  there  is  no  charge  per  length 
of  weld.   The  installation  cost  for  the  shear  connectors  is 
50  cents . 

Design  the  girder  until  the  change  in  the  critical 
moment  is  less  than  one-half  percent.   Determine  the  deflec- 
tions at  the  centerline  of  each  span  and  use  the  debugging 
output  option. 

The  input  data  for  this  problem  are  shown  in  Table  C-l. 
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Problem  No .  2 

Design  a  noncomposite  girder  for  a  70  ft  simple-span 
bridge.   Design  the  girder  in  A36  steel  and  for  an  1120  AASHO 
loading.   Use  a  load  factor  of  75  percent  and  do  not  include 
impact  in  the  design.   Read  in  the  coordinates  of  fifteen 
analysis  points. 

Optimize  the  girder  design  based  on  weight  and  use  a 
total  of  three  design  cycles.   Determine  the  deflections  at 
all  points  and  print  the  solution. 

The  input  data  for  this  problem  are  shown  in  Table  C-2. 
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001257  D'1  -  COuSP(J,3,T|  -  rONSP(Ji?iI) 

O01?h*  NSP  =  OM  /  C0NbP(JiI,I)   ♦  ^AF 

00127?  SP  J  DM  /  NS>P 

001?74  UP  1  =  SP  /  ■>. 

001277  SP  =  MP3  »  3, 

001301  C0«J"»PCJi  1*1)  =  SP 

O01303  NT  a  NT  »  Dm/SP  -  MSP  ♦  NAF 

00131?  If  (hT  .RE.  NAfni  (50  TO  40ho 

001314  J  =  J  -  1 

00131*  IF  u  ,F(],  JK)  r,n  TO  40on 

O01317  00  IC  151 

001320  iso  n  =  (.()nsp(  jk.t,  1 1  -  ronj(j4) 

O013P5  NSP  =  I)  /  CUNSPl  |K.l  ,T)   .  NAF 

001333  SP  »  0  /  NSP 

001335  NP *  i  SP  /  3. 

001340  SP  *  MP  1  «  3, 

00134?  CON  >P  (  JK«  1  >  I  .  II   r  SP 

00134*  CONiP(  JK«1  ,3,  J)  -  i-T-iSP  I  |K  .1,  I  I 

00135?  CO^ih  (JK«  I  «?,  I  )  =  ")NSP|  ik»),3iI>  -  "ISO  •  So 

00135*  COkIsF ( JK,3iI)  =  r.OMSPI  JK*-I«?.D 

001 3*1  4="0  CONI I MUF 

001  3M  40"n    CON  I  1  Nl  t(- 
0013*4  0"    «P     I     r     1,     >■■>■ 

0013*5  KSC1      =       I 

0013**  00     .1      )    =     1 ,     i < 

001370  IF  (ARS  (rt)NSP(  J,  1  i  II  1  ,|T.  .oOOOl)  RO  Tn  30 

00137*  n  =  CONSP <J. 3. I)  -  CnNSO(J,?,I) 

001403  NXO  =  D  /  C"'J<  I  |  |.  l  i  1 1 


CONO 

'4? 

CONO 

?»3 

CONO 

'44 

CONO 

?45 

C0Nn 

?4fc 

CONO 

?47 

CONO 

?48 

CONO 

?49 

CONO 

?S0 

CUM) 

?S1 

CONO 

?52 

CONO 

?S3 

CONO 

?54 

C0Nn 

?55 

C"Nn 

?So 

CONO 

?n7 

CONn 

'5H 

CONO 

?5<> 

COM) 

?*0 

CuNO 

?ftl 

CONO 

?*? 

CONn 

?ft3 

CONO 

?fc4 

conn 

?65 

COM1 

?ftft 

CONO 

?*7 

C  'NO 

?*fl 

CONO 

?f)9 

CONO 

?70 

COM, 

'71 

CONO 

?7? 

CONO 

?73 

CONn 

?74 

CONn 

?75 

COno 

?7ft 

CM) 

?77 

CONn 

?7« 

CONn 

?79 

CONO 

'«0 

COM) 

?«1 

CONO 

?B? 

CnNO 

'P3 

CONn 

'R4 

CONn 

3P.S 

CONn 

?fl* 

CONO 

3B7 

CONO 

?«S 

CONn 

?fl9 

CONO 

5<90 

CONn 

"»1 

CONO 

?9? 

CONn 

?"J3 

CONn 

7Q4 

CONn 

pq<? 

CONO 

?"»* 

CONO 

?<>7 

COno 

Tin 

CONn 

?qq 

144 


p  0 1  *Oft 

If     I 

ooi *  i  s 

If     ( 

00l*?l 

KSCI 

ool*?? 

IF     I 

noi*?* 

If     I 

001*30 

COM-. 

noun 

If     I 

ooi*«o 

Gn    i 

001**1 

4? 

Com-i 

001*50 

COM-> 

001*51 

GO     f 

001*51 

43 

COMj 

001*63 

COm» 

001*66 

*1 

COM! 

001*70 

39 

IF     ( 

001*71 

on   •« 

001*73 

D1    '» 

001*74 

CI  It 

001504 

4S 

COM! 

0015)? 

KH    * 

001513 

DO    •» 

00151S 

IF     i 

0015?3 

IF      ( 

noi5?4 

IF     1 

00153S 

CO-MI 

0015*1 

(SO    f 

001541 

'.7 

com-. 

001547 

com-. 

00155? 

CO-M-. 

001555 

KH      = 

001556 

4  6 

CO'MI 

001560 

40 

r 

COMI 

001563 

DO     '» 

001564 

DO     4 

001565 

IF      ( 

001573 

xn    • 

001600 

NXi) 

001603 

IF      ( 

00161? 

xn    - 

00161h 

If    < 

0016?? 

CO  *  ■> 

0016PS 

COM-. 

0016?7 

Git    I 

0016?7 

400? 

If  l  • 

00163S 

COM-. 

001 644 

CO  M  i 

001647 

GO     1 

001647 

40'<3 

CO  Mt 

no  1  6^ 

*V'» 

CllMI 

001663 

If       1 

001665 

no  » 

00166h 

j  i  = 

001 667 

if    i 

00167S 

4I1US 

CO-MI 

O0167f- 

4nofc 

J  )   = 

001 700 

XI)     - 

001 704 

DO    ■• 

AHS  (NXDOCONSD  (  1,1,1)   -  n>   ,GT.   .00011  Mvn  I  MXn  ♦  1 
IjXD  ,GT.  3)  RO  TO  41 

=  ) 
J  .FO.   I)  1-0  10  4? 
AHSICONSPIJtl  >  1  ,  T  I  I   .LT.   .0001)   GO  Tn  41 

P<Ji  1 .  I )    ■    rniMSo(  j-1,1,1) 

CONSPI  J«l  ,  1  ,  1  >  ,LT.Cn>MSP<   1-1  .  1  .  T  )  )     fn-isp,   i,  i  .  1  )  rCONSP  i 

0    41 

P(1,3.I>     =     rON5P(l,?,I)      «     1,     »     (-ONSPM  ,1  .  Tl 

P  I?,?,  I >      c     fOMS?  (1,1,1) 

0  41 

P(J.?.l)     =    rohSJ ( j,i, I )     -    3,    o    rONSP (  1, 1  .  T> 
PlJ-I  ,3,1)     =    COmsPIJ,?.  i) 
J  NUE 

IvSCl     ,E0.     0)     (id    TO    40 
5    .1    =     1  .     1 6 

5  K    =    It    3 
h , I  I    =    CONSPIJtK, 1 ) 
F (j , k , I )    ■   n . o 

l 

6  J    =     1  ,     1  6 

AHSICIJtl  ,1)  )      ,LT.     .0001)     GO     TO    u() 
J     ,1').      1  )     GO     TO     l,  I 

ARSic{j-i,i.T>    -  ci. it )«!)>    ,gt.    .onnTi    r-o  Tn  47 

F  («»--l  ,3,  I  )     =    r  I  1,3,1) 

G    4/. 

f  («",  ltii   =  r  1  1,1,11 

f  (K  t,?t  1 1     *    r(Jt?t  I) 

P(k ■■!,  't  I  1    =    r  I  ).i.  I  ) 

Kls    »     1 
IN  US 
lNIII 

INAI   AOJIJSTMFNTS  OF  THE  fONNFCTOP  SPin'r. 
300  I  =  1 ,  MM 
39n  .1=1,16 
AHS ICONSP ( J, 1 , T ) )  .LT.  .0001)  Go  To  430" 

CO.ISPI  J,3t  I)  -  CONSP(Jt?tI) 
=;  X'   /  Cic  SP  (  l.l.p 
AMS  IMXO«CONSp  (  I,  |  ,  I)  -  XOI   .GT,  .000)1  wvn  5  NXO  ♦  1 

MXn  *  CON'SP  (  J,  1  ,  I  ) 
CONbP  (J,  1  ,  1  ,  1  )   .|T.  POMSP  I  J,  1  ,  1  )  1   GO  Tn  ■  <iq-> 

p(j.3,D  =  rnN53(j,?,n  ,  xn 

P(  J»l  ,2tl)  =  CONSfl  I,  3.  I ) 

1  4  IPO 

eS(l   i.nSh  (  J*  I  ,  l  ,  l  )  )      ,[T.      ,0001)      GO     Tn     '.or-, 
I' I  1,3,1)     *    rr)NSO(  l,?,l)     ♦    Xo    -    rOMSP  1  i.i.Tl 
P ( J* 1  ,?« 1 1     =    ro\j«^p  (   1,3,1 ) 
(     4  •'"' 

p  1  1. 3 , 1 )    ~   ro.vSJ  1  i,?,  11    ♦    xo 

IN  Iff 

I'M      .  f "'),       1  )       RO      TO      1  1 

IJUS     I    =    1  ,    1 1 
J 
AHS  t  r   iimsp  i  .1,  -.,  mvii  )     ,|_T.    .ooni)    on    Tn    4nu< 
INUI 

=   JJ  -    I 

I  F_   .:,  I  1     -     r  0»  HP  (    I  ),  1,  m  II 

•  I  9  r   .1   =    I  ,     II 


CONI) 

300 

cono 

301 

rofMD 

30? 

COM) 

303 

COM) 

30* 

CONI) 

30S 

1*1 1 

1,1)     CONO 

306 

cond 

307 

CONO 

30R 

COM) 

304 

f  ONI) 

310 

CONI) 

311 

CONO 

31? 

CONO 

313 

CONO 

314 

COMI 

315 

CONI) 

316 

CONI) 

317 

CONI) 

31A 

CONO 

319 

CONO 

3?0 

CONO 

3?l 

CONO 

3?2 

COM) 

3?3 

CONI) 

3?4 

CONIl 

3?5 

CONI) 

3?6 

CONO 

»?7 

CONO 

1?« 

CONO 

3?9 

CONIl 

330 

CONO 

331 

cono 

33? 

COM) 

333 

CONO 

33* 

CONI) 

335 

CONO 

136 

conii 

337 

CONO 

33H 

COM) 

339 

CONO 

340 

CONO 

141 

rONn 

14? 

CONO 

3  4  1 

CONO 

144 

CONO 

345 

CONO 

146 

CONO 

34  7 

CONO 

1<,H 

CONO 

149 

CONO 

ito 

CONO 

151 

CONO 

1^? 

Cono 

"<^3 

CONIl 

>S4 

C   INO 

IS^. 

COM1 

1=.h 

CONO 

1G7 

1  ].', 


00170ft 

CO-IlP  < J. 3.MN) 

=    r'i'iSp  (  i 

, I.mni     ♦    »n 

nomi 

40R7    CO  M->R  (  J.?.MN> 

=     <-<>'lSP(.J 

,Pi»N|     ♦    Xn 

ooi7?o 

1  n     CONI  lNUF 

0017?0 

RF  ruHN 

ooi7?i 

END 

COM)  3SH 

CONO  35<J 

COND  IftO 

COM)  161 

COM)  3ft? 


SUHPRO'iRA^I     LfN'JTM 
002531 

UNUSED    (-U--1RII  FK    SRjCF 
034100 


1 
? 

00000?     "                     CO-MUM    ,JAHE  (  1  8)  .SARFa  (  1  h)  ,DfpTm  (  1  B)  ,Fl   <in  (  l  cm  ,Fi  TH  (  1  R)  .wthTmI  18)  i  CPUS  3 

4 

ft 

7 

^f  iniwr  ->nioii*f^i  fiics'iijM  i  ri  i  ibi  nir  ^"<-  a  (  ^  ,  4  >  ,  "fssh i  iftvf  *.  j  ,iraHi^tji  (Lkii^  8 

ftCOM  (R)  .Consp  (  lft,  3.4|  ,Cn\yPi_  I  1  ?»ft.?l  CPns  R 

nooop?                    co i iON/une /jpT3, ipt i i  CPns  in 

00000?                          COMmON/TwO/lPCM  CPnS  11 

00000?                          CO-MON/SIX/FU  CPOS  1? 

00000?                        CO»iil/N/t- H.HT/TST  CPnS  13 

00000?                          CO<Ml>N/(EN/HAUNCH  C^DS  1* 

00000?                          COM"tON/ElEVEN/Kra,KCt«;iS»SHStSLS»I5S,  I.  SS.HMn  .JITST  CPOS  IS 

ooooo?                  clmiOn/  ifen/pth.kt .  i  CPns  ift 

00000?          KF4I.  I*i  It  .JSIFFi  ,  rLCON,  IMCON.MR.LFNfiTH  CPnS  17 

00000?          DIlFNSIDN  LPC"  (  1  0)  .PTH  (40  .4)  ,RSm  (40  ,  4  )  ,phm  ( ,,  ~  ,41  CPnS  18 

00000?          IF  (1PT10  .OF.  3)  write  (h.lo4<))  CPOS  \9 

ooooii    loco  format  (  imi  ,?ox.  iihstart  cpofS//)  CPns  ?o 

C         ^fcRO  OUT  Tl-F  RF'JJTRFn  ARRAYS  CPDS  ?1 

000011          DO  w400  1=1.1?  CRns  ?? 

000013  no  v4oo  j  =  l,  ft  CPns  ?3 

000014  do  <4oo  k  =  i.  ?  CPns  ?4 

000015  R400  COVPL  UiJ.K)  =  n.O  CPOS  ?s 

000031                    mn   =  ns  -   l  CPns  ?ft 

000033  J    =    2  CPnS  ?7 

000034  ft30(i    DO     <(>01     I     =     1.     mf  CPDS  ?8 

oooo3ft        no  ho)  jj  =  i,  4  CPns  ?r 

000037    9oM  pr-iii.jj)  =  o.o  CPns  30 

000047                          DO    M'ftl     I     =    1,     40  CPDS  31 

000051                    no   -.eft)    k   =    1 ,    t,  CPns  3? 

00005?           5(ih)    pmil.M     =    p. 11  CPns  33 

ooooft?                   km   =    1  CPns  34 

0000ft3          NC  *  KH  CPns  35 

r           DETERMINE   IhF  i'H  ARRAY  CPDS  3h 

ooooft4                   no   -iiono    I   =   1  (   nf  CPns  37 

noooftS                       IF    (I    ,1.1.     11    r.o    1  1    rjson  CRns  3H 

000070              RS10     Pf-H(K8,1)      =     P|   Mi'll,    II  CPOS  3R 

00007S                            PT>1IM,J|     =    PLwTDd.    I)  CPOS  40 

oooioo         PTHiM3ti)  =  cnoRii)  CPns  4i 


1  If. 


100101  PTHOHitl  =  C0OB(I«l)  cpos  4? 

000103  GO  IU  9i  in  CPDS  43 

000103  95nn  IF  I  <>H  S  (  PL  A  T  H  I  I  ,  I )  -P  TH  (  KR  ,  1)  )  .  L  T  .  ,  f|  0  0  0  1  I  On  .fl  0600  CPDS  44 

000113  kh  »  KH  ♦  1  CPnS  45 

000114  NC  -  KH  CPDS  46 

000115  SO  Id  9511  CPDS  47 
000115  9610  PTr|(K.H»4)  =  COO«(I  +  l)  CPnS  4fl 
0001?0  9000  CO'<l I  lNlifi  CPnS  49 
O001?3  KSS1  =  n  •  CPOS  so 
nooi?4  IF  inn  «F.Q«  ])  r.ri  TO  7no'  CPns  51 
noni?s  if  (j  .i.ij.  li  on  in  700?  CPns  5? 
oooi?7  no  'ooi  11  =  ?,   us  CPns  53 

000131  1  »  10SHI1I)  CPOS  54 

000133  DO  /(.0  1  JJ  =  1.  NC  CPOS  55 

000134  IF  ICOOH  I  I )  .fiT.PTH(JJi31  •ANO.COORI  I)  .LT.oTHl  lJt4)  )  r,n'TO  7005  CPOS  56 
000146  GO  fO  7103  CPns  57 
00014S  7015  10  *  JJ  .  CPDS  5fl 
100150  IF  it'TMII'J-1  ,1  )  ,LT.  .1)  GO  TO  7006                                   CPns  59 

000153  KSSI  =  I  CPns  60 

000154  PTHIlO-l,l|  s  0.0  CPnS  61 
0011155  PT-UlO-li?)  =  r.n  CPOS  6? 
000156  7006  IF  IHnI  l'J«)  »  1  )  .I.T.  .1)  GO  TO  7003  CPOS  63 
000161  KSSI  =  1  CPOS  64 

00016?  PTriiln«i»i)  =  o.i  CPns  65 

000163  PTH(IO*1«?)  »  fl.n  CPns  66 

000164  70^3  Cf>"tfINUf  CPDS  67 
000167  70H  COMflNUI  CPOS  6R 
000171  IF  (KSSr  ,E().  0)  GO  TO  700?  CPns  69 
00017?  Kfl  =  NC  CPns  70 

oooi74  ktt  =  i  CPns  71 

000175  NC  =  1  CPOS  7? 

000176  no  mho  i  =  i ,  kr  CPns  73 

000177  IF  (I  .i.T.  1)  CO  TO  7015  CPDS  74 
000?0?  70?0  DO  /01?  KKK  r  1,  4  CPDS  75 
000?04  701?  PT-MKTT.KKK)   =  DTHIt.KKK)  CPDS  76 

000?15  go  ro  7nl0           /  CPns  77 

000?is  7015  IF  (ARS(PTH(I«1)  -  »THU-1«1I>  .LT,  .oooil  m  Tn  701*  CPnS  7R 

noo???  k  t  r  =  kit  ♦  i  CPns  79 

ooo??4  nc  =  ktt  CPns  so 

000??5  GO  ID  70?1  CPns  fll 

O00??5  7016  PTniKTT.4)  =  PThiI.41  CPns  fi? 

O00?30  7010  CO^MINUr  CPns  A3 

r  ->k!Nl  THK  ri-ir.l'-Ai  ptm  uroAY  if  dfootoFH  Cpns  B4 

00OP33  701?  IF  IJPTln  .LT.  3 )  RO  TO  )oo  CPns  A5 

noo?36  wPiit  (6.ioooi  •  CPns  as 

000?41  lonp  FO-l-iAT   (  1  M0>  10»  ,  1  HH-IDIGIHAL  PTH  ABO/tV//)                                  CPnS  A7 

ooo?4i  no  iooi  l  *  i .  nc  CPns  rr 

O00?43  mil  WW  f  I  t  I",  100?)   (PTHtT.K),  K  c  1.  4|                                      cpns  R9 

000?61  ion?  FO-MAT  (?Xt1HF'l»TE  Is  ,F10.5,3H  X  ,F6.  f,  KX  .nHS'TART  e  ,Fffl.?»lflX«  CPns  90 

I6mf  <o   =    .f in.?i  CPns  91 

noo?6i  inn   k'>   a   n  CPns  9? 

000?6?  KILL  =  :'  CPns  93 

000?61  KPST  =  0  CPnS  94 

000?64  IF  (NC  ,FQ.  1  .AMD.  AMS(OTM(),H)  jt,  .nnni,  r,o  To  10?  CPOS  OS 

C  -,1AUI  THF.  lono  rip  ihc  oLftTE  CUTOFF  CPns  96 

000?76  DO  J300  1  =  1 .  MC  CPOS  97 

O00?77  KST  =  Kl<  CPOS  9R 

00030)  KB  =  n  CPOS  99 


L47 


00  0  30 1  *T     »    li 

00030?  IF     (IPIir     .,F.     31     -  ^*  T  TF     (  ^. .  1  <">  *►  1  >     I 

00  0  3  1  1  ln'1     F0rt«Al     (  /L*X  .  ?  ]  pr  n  I  r  l|   >  T  T  O'jS    FOk    I     n     .1?) 

000313  IF     iF'TMIItll     .LI,     .1)     SO    TO    0)00 

00031*,  if    it  li  l    .f  i.    ni    fin   to    i9o 

000317  IF     (AHSIftHt  1 ,1  )-PTH(J-l  1 1)  1     »GT.     .oOOl)     fifl    TO     '90 

P003?S  30     Hi    9300 

O003?S  1P0     IF  (  >ns  ll'THI  I  .  1  1  -pTh  i  J-]  .  ))  )     ,l.T.     .orinl)     «0    in    9">o0 

000335  IF     <i\sT     ,KU.     ?)     RO    To    9">S9 

r  itTFi'MNt    1HF    PlATF    FND    CONDITIONS 

O00334  DO     »33d    K    =     i  ,     bja 

00033*,  M    =    K 

000337  IF      (  A.MS  (COOW  (  I  A) -PTH  (  I  .  1)  )  .1  T.  .(10011  I      RO     Tn     "til 

000 34S  9330    CONI  INI1I 

000347  933)     K'J     -     IA 

0003S\  KC     =     1 

0003S?  IF      (I     .i    j.     1  )     Pn    TO    9  133 

000 3S 4  IF     (PlMll-1,1)     .OT,     .1)     Kf    =    ? 

0003157  9333     ISS    =     1 

000360  IF     (Plll(I-lil)     ,r,T,    PToiT.lll     ISS    =    ■> 

000364  IF      (IPIIF      ,£().     4)      " « T  Tt-       (6.101?)      K'l,     KF,      ISS 

000*00  101?    F(M«AT     i/S<,k^ii    =    .i?,?X,ShkC    =    ,  13, ?x  ,<,Htsc;    a     ,I?| 

000400  IF      IMN     .Fu,      |)      fin     TO     ^3?9 

r  '(ETfc'iiMlNd     Tmf    maximum    ALLOkAmI  F    rilTOFF 

00040?  DO    «310    K    =    ?.    mm 

000404  10    '     lOSP(K) 

000406  ]F     irP-,1     .£Q.     II     RO    TO    7701 

000410  IF     (I  no.j  (  ID)  ,(F  .pTh  (  j  ,  3)     ,Arjn.    COOf>  ( 1 0)  .  I.F  .  P  TH  U  1 4 )  )     GO    to    93|S 

OOO1"??  GO     10    4  <1  0 

0004??  77"0     IF     ((Oil.   (ID)     ,r-F,    PTh(I-1i3)     .AIiP.    roimiiH)        IF.     PTH(I-1,4)) 

1GO     fO    4317 

000434  931  0    CP«J  I  I  Mill 

000437  93?9     IF     (6  1  LI      .El).     1)     GO    TO    03]7 

00044  1  DMA*    =    SUHLfN  (Ml 

009441  IF      (ISS     ,FU.     ?)      ru-lAX     a     S'lRl  F  N  I  I  A- 1  ) 

000447  GO      It)     9(?1 

0004SO  931s  IF  (KH  .  Fi.  ?|  r,o  TO  9317 

0004S?  ND"  =  OFPTH(IM)  ♦  19. 

0004S6  0  =  Mil)  »  1.0 
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ooooiss 
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oooi  i  ? 

»  1  i  >l  = 
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fS  FOP  Thf  IMiTial  VALUES  IN  ThF  Oil  -/AT-JX  VhIcH 
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?ori 
?9r4 


SUi-<iiiJ 
i  if  T 
COM  10N 
Dl  Irr.S 
IF  IIP 
F  O-ciAT 
DO  1^0 
IF  ill 
PL„  = 
GO  10 
CC N  I  IN 

00  <70 

SITrld 

iiET 

si  r-.(  i 

K  a  C 
IB  =  1 
DO  v7  1 
DO  47? 
IF  I  AH 
CONI IN 

si  t-kk 

K  -     K 

19  s  k 

CONI IN 

HLA 
CALL  S 
IF  I  IP 
<K  |  I  t 

FO-eiAI 
DO    -<M 
Wrt[  IE 
FO-l-iAT 

no   ^o 

DO    -"*0 

rpi  1 1 j 

DO    ^"»0 
ill  T 

pcr«i<i 

1  =   I 

IF      I  AH 

I  =  1 


I  HE    Fl.OPT     (PTn.Nr.COST.FT) 

IMTZE    THF    iotatIO'iS    OF    ThF    Cm/FR    Pi^Tr    sp,  iCFS 
/uNE/ JPT  ?, [P r )  i 

ION    STTHUn)  , ?TH  I  4  0.  *)  , PC  TS ( *n )  , rnsT (o.  ,ri|«(ii  ,TP("1,*0) 
TIO    ,r,p.    i,     «9tTE(AilOP) 
<  1H1 1  in*  ,  ivhFnTky    TO    FLUPT//) 

0  I  =  1  »  HC 

Hiii?)  ,lt.  *.o>  so  Tn  ivoo 

PfHUt?) 

1V0? 

UF 
[    >    li    KC 

1  =    0.0 

fHMlNE  THF  STTh  ARRAY 

)    =    PTH(I,i) 


I    =    2i    NC 

KK     =     1 ,      TH 
SlSTTH(KK)     -    BTHUiUI  .LT,     .00011     f,f    T"    971 

UE 
)     =    PTHI 1,1) 

♦  1 

-     1 
UF 

CF   THE  STTH  ARRAY  IN  ASCEMOINfi  DRDP^ 
OkT  (STTh.TH) 
THi  .LT.  3)   (50  TO  400 
(6, 4H0) 

( 1  ho, I o* > 1 5HSTTH  aRRay //) 

1=1,1* 

<'.,9b?)   I,  STTnlII 
<  1  OX, I?, si , Fin,*  I 

0  1  =  1  .  Mi" 

0  J*  =  1  ,   TH 

K)  =  o.r. 
11=1.1° 

FriMINE    THF    STARTING    SECTION 

1  =0.0 

s  (PTH<  i,i)i         ,r,T.    ,nnni)    fin    Tn    ?o^3 

♦  1 
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ooou? 

000U7 

000151 
00015? 
000160 
000161 

0001M 

000161 
000165 
000170 
000313 


00021  3 
000215 
000220 

000222 
00022* 
000226 
000232 
000236 
000237 

0002*0 
0002*2 
0002*7 
000253 
00025* 

000256 
00026? 
000273 
0  0  0  3  0  1 
000306 
000315 

000315 
00031 7 
0003?* 
000331 
OOO  3*  1 
0  0  U  3  *  ? 

0003*3 

000350 
'000356 

000  156 
0  0  0  3  6  ? 
000  3f>h 

o  o  o  i  n 


no   10  ?90* 

290J  lb  ■=  I 
C  0tTE><*IN£  THE  FNUlNU  SECTION 

I  ■  NC 
290b  IF  (AFtSlPTHlI.il)     ,01.  .0001)  1,0  TO  2905 

I  •  1  -  1 
(30  TO  HOh 

2405    IF    3     I 
C  HE6IN    TIE    OPTIMIZATION    OF     IMF    COi/tk    PLATES 

00    *0OO     I     ■     U.     \f 
PLl EN    ■    PTHI 1 ,*)     -    PTM( l,3i 

IF     lIPTIv     .£0.     *)     wollt     (6.101)     IS. it. I iPLW.PLLFN 
101    FOHMAT  1  1H0,5H!S    =     ,I2»*»5hIt     •    .I2»A.*MI     ■     ,I2i*i6HPL*    *     »F5.1i 
IX.HriPLLtN    =     ,  M  0.?) 
C  INITIALIZE     In!     STAFtTlNG    CONSTANlb    AND    COSTS 

1)0    *0Ul     K     «     1  ,     1H 
*001     cl  IK)     *    0.0 

IF     (  1     .£0    .     IF)     r,n    TO    4500 

c       vapt  rut  Thickness  of  thf  ao.jacen!  section 

DO  *01 0  *  «  1 .  1H 

TCT  *  9  '999. 

IE  (STTrt(K)  ,r.T.  PTm  (  I  *  1  *  I  )  i    GO  To  *020 

IE   (AFiSISTTH(k)  -  PTH(1*1»1>)     .LT,  .000011  GO  TO  *0?n 

CT  (M  *  0.0 

00  10  *0  1 0 

c       vaht  me  thickness  in  thf  oesion  section 

*02u    00    <,030    KK    =    1.     in 

IF     ( ST  Til  (KM     ,1,1.    PTH(llll)  GO    TO    *0*0 

IF     (AHS  ISTIrtlKM     -    PTHtl,!))  ,Ll,     .00001)     GO    TO    »n*0 

(,0     10    *U30 
•  0*0     IF     (1     ,nf.     IS)     t.O    TO    *0S0 

C  OETtrtHlNE     THE     COS!     OF     THE     INITIAL    SECUON 

Tma*    *     TOPI     (SlTHlKKl ) 

CEk    e    CU5T  (b)  *H  w»fOST  (U)  ♦  ,5«TMAA»«<:,PLW1'CnsT  (  7) 
rbw    =    2. "COS)  U .) .?.•       HLLEN       "(COST (H) «.01759*»C0ST (71  ) 
IF      (AhSlSTTniK))     .CT.     . OOOOUll     GO    Tu    *0*1 
Cbw    -    C^'ST  (bl  .Pl.w»COSl  (U)  »PLW"TMAA»»^«COSI  (  7) /2. 
CO     TO    *>l*2 

C  lit  H'lilNt     thF     msi    UE     THF     BOll    Wt-Ll)    IE    heooi^eo 

*0*1      1H     •     SI   IHUI 

IE      (SIT"  (KM     ,LT,     SllulK))     Th    ■    STTlllKKl 

]E     (AE)S(STTri  (KM     -    STTh(Kl)     .LT.     .Ouudl)     Go    TO    40*1 

CHw    =    COST  lb)  .  Flw«  (COS  MR)  ♦  ( ,  12S*TH»,2H8676i»TH«Z»  »COSTI  T) ) 

i,()     10    *  )*? 
*0*3    Cm    =    D.O 
C  nlll  '<-U  \it     Tut     maTEH1«L    roSl    Of     TuE    Pi  Ale 

*0*2    C'T     =  PLLEN       «Pl.w»STlh(KK)  »C0SU1)  «*90./172rt, 

C  liETfHMINE     |n|      KHAL    COST     UE     1HF     btCTlON 

re     -    Cel     .     Chw     .    CF  ti     ♦    Lsw     .    PCTSIKM 

go   ro  ti/i 

C  Hi    lC'«^lNt      Thl      rnsT     OF      Ti^      bF  AL     HtLOS 

*0bu     IF      (Alib(STTrl(Kh)l      ,|   I  .     .  000  00 1  I     i>0     10    .,0  55 

II  lAMbtsl  iriir  I  I     ,i  1,    .uuonii    GO    10    *0S6 
CSw    *    .'.*       PlLfN       »  ICOSl  (Ft)  ♦.017b94»C0ST  (  7)  I 
r ' «    ~    0.0 
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000374  t,0     10    4  i4  1 

000374  4055     If      lAbSISI  lrt(K)  1     .IT.     .U00O1)     tin     Hi    »u'j7 

C  UF1fcXM|NF     In)     fOSl    OF     ThF     APPmOUmATE    FND    WFLOS 

000400  J»H    «■     I  OP  I  (SIlH(K)  ) 

000403  Ct*    =    COM  (o)  «3.  .Pi  w»tUbT  (H)  .FLW".5°  I  max  ••-"'COS  I  17) 

0004 ib  rc  a  cf «  «  pctmkk) 

00  04?ll  (H)  TO  4  J  70 

0004?)  405b  Tmia  =  IPPI  ISTlli(KK)  I 

0004?5  CtS  -    C  JS  I  (  M  .n  w<.f  uS  I  (o)  ♦  ,s»  I  MAAu"2BCOST  (  7)  »PL« 

00U41*  C^w  *  0.(1 

000417  rSW    =    r".°       Pi  I  f.  N       «  If  ob  I  (ii)  «  .0  l  7S^4«COST  I  I )  ) 

0004  4  3  (iO    TO    *u4  2 

000444  40">7      FL     =     I'CISHM 

00044*  40'0     If      (IPI1II     .to.     4)     *-Mlk      lb. 102)     STTnliv),     STTh(Kk).     Tr 

0004*5  HlZ    FOwiAl      O*  .  t'lHl  ill  TknFSS    ON    KlLnT     -     iE  lO.b |2A i20HTHlf KNFSS    ON    LEFT 

1  =     ,F  10.  >»?»i  'H(  (1ST     =     .Fill.?) 

0004FS  U      111      .<",!.     Irll     r,o    10    4030 

C  SlOi't     Trir.    llisl     r,h|)     UiF     T>llLI*.NtSS 

000471  Tl-I     =     It 

000471  |M'    «    M 

00O473  40  30    fONl |NUt 

C  SIOHf     lot     THICKNESS    KFODIHtU     IN     l^f     IP    AHHAY 

00047b  TP  (  I  »K  J     =    b  I  in  (  it-Pi 

00050?  CT  Ir  )     =     TCI 

000504  It      (IF-TI"     .hi,     4)     VMlIt      lb. 1101     UUI  i     TP(l.K) 

0005??  110    foi)maT      (  IrlOt  lnX.7HC')Sl     =     ,F  1  b.b.bAtbhTP    =     .MS. 5//) 

O0U5??  401  0    C^'j  1  IMIE. 

C  DE  TFi<m|  s£    THl     Pfls    a*kay 

ooosps  no   40^n   **<   =    i  .    |m 

00OS?h  40*0  PCTsmuM  *  CI  ^M(l 

000532  IF   IlKllO  ,UT.  41   l,(l  10  40011 

000534  Wh  II)   C>,  4  1  0  I 

O0OS40  4lU     FOllMAl      I  lMll»  lflhPCTS     AHHMY//I 

000540  wHTTt      CmIOJI      (f'CTS<KIMOi     k*K    ~     |i     ID) 

0005^3  10J    FOI'hA  I     I  4  I  M  5.5  .'■-X  )  ) 

0005S3  4000    COM  lNUt. 

C  vAUt     Int     IhlCKNFSS    0)     1  HI     LASI     StCTION 

0005*0  4500     00     tlUII     K     s     1,      (H 

0005b?  U      ISTTnln]     .il.     PlHIIt.ll))          00     10    »l?n 

0i)05h7  If      l»HbtSrtrt(K)     -    Plullt.ll)           .LI.     .000001)     00    TO    4I?0 

000573  PCTSir)     =    0  .0 

00  05  75  GO    T.I    *  11! 0 

C  DF  TF.x'MNt     T  mF     rnsT     Of     ThF    LASI     StCTION 

000575  *U0    Cff     =    SllHlMo          Pi  *             o    PILtN       "          CliSl  i  1  )  o*iJi)./17?H, 

000*0?  C^w    =    ?«"       PlMN          "  (COST  (HI  »  .01  7b9,<•(.0<;,  I  7)  I 

000*07  If     (If     .f_i.     I-)     CSV.    «    Lbw    .    ?.    •    COjIibi 

000*1  4  T^x    =     I  0P1  I'M  lh(K|  I 

000b?0  CI*    c    COST(t>).             P|  «             »r(»Sl  (H  )  ♦  ,'j"  IMAX»«?«COST  (  7  |  o       Pi  w 

000*31  PCTSli-l     =    BCTs(K)     .    CPL    .    CSW    ♦    Ltw 

000*3*  ]F     lit     .ED.     ISI     PrTS(R)     =    PfTSlKl'.    Cf  « 

000*4?  4ioo   continue 

000*45  rCT      «     9'<99^  . 

C  lip  1  F»<M  J  ME      iHfc     THICKNtSS    of     ThF     L«bl     SECTION 

ono*4h  no   4 1  in   «    =    i  ,    in 

000*so  If   (AHSlPCISlKl  ]   .IT.  .ooonnil  r,u  lu  Ml" 

000*54  IF   IPC'^I")  .1.1.  T  r  1  )  GO  TO  4130 

000657  TC  T  =  P(.TS(<) 

000**0  ICt'  c  k 
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FOPT 
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FOPT 
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FOPT 

130 

FOPT 
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FOPT 
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FOPT 
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FOPT 
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FOPT 
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FOPT 
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FOPT 
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FOPT 
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FOPT 
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FOPT 
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FOPT 
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FOPT 
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OOOftftO 
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0006ft") 

oooft6s 

000hft7 

0U0ft73 

4bO 

000h73 

000700 

5520 

0007?) 

ssai 

0007?! 

0007?7 

435 

0007?7 

00U733 

5522 

000751 

552  3 

000751 

420 

O007M 

0007S* 

0007S5 

001000 

<tbo 

c 
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00100? 

00 1001 

41  10 

001005 

00 100ft 

001013 
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S|)OU 
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41SU 

0010  3?. 

00 1034 

o  o  1  n  <*  ? 

no  lius 

00  104ft 

0OlO4ft 

4  1  bO 

0  0  1  0  <i  h 

o  o  l  o  s  i 

oo 1 ns* 

linu 

0  0  10  54 

oo loftn 

1 11' I 

ooi  1 14 

1  1U2 

001114 

0  0  1117 

1  »4 

001117 

4^.1 

n  u  1 1 1 7 

001 1?0 

CONTINUE 

IF     (1PT10     .LT.     3)     GO    TO    420 

11     a     It     -     1 

rfhlTf     (h,4t>0) 

FORMAT     HHOtlOXi JBHFINAL    TP    »HH*Vj...//| 

1)U    5520     I     *     Is,     IT 

WRITE     (6t5bil)      I.     (     rHll.M,     K    ■     I«     lb) 

FORMAT     (5X,l2,5x,10(F7.*,3X)  ) 

WHITf      (ft«455) 

FORMAT     UH0ilnX,3HHCALCUl.ATlONS    FOR     IHE.LAST    SFG*FNT //> 

riO  552?  *   »   i ,    ir 

WRITF     (ft»55iJ3l     K.     S  I  TM  tlV  >  •    PCTSIM 

FORMAT     15X,*hi\    =     , |?,2X» 7HSTTH    o     iF1Q,5,2A,7hCOST    ■    ,FI5.5) 
CONTlNUr. 

■    Plh(  IKiI)     =    S1THUCHI 
Fl     «    Pcrs(lCH) 

IF     (IPTin    ,GF.     3)     WHITE     (6,4h0)     It.    STTHUCHli    FT 
FORMAT     ( 1HU »?1 MlHlTKNtbb     In    SEGMENT     i12,JH    o     ,F10.5,|5h    WITH 
IT    »     »F  1  >.">) 

DFTtR«|Nt    THE     THlCKNtSSFS    OF     TMF    HEMAlNlNG    SECTIONS 
I     =     If 

IF     (IS     .Fi).     IF)     GO    TO    *lftO 
l)U     4  140     K     t     1  ,      |H 
ICR    -    K 

IF      (M*SlSrirt(K)     -    PTH(lil))  ,(T.     .000011     GO    TO    Also 

COM  I  MJt: 
wR  1  1 1      (",5000  ) 
FORMAT     <  1  HOt  1  OX  i?0MF  HHOH    FOUNU    IN    FLOP!) 

slop 

PlHI  1-1  .  1  )     a    TP(  1-1  ,  ICK) 

H  Tm ( I  -  1 »  ? 1     »    PL  * 

IF      (APVfTrtM-ltllJ     ,LT.     .00001)     PTn(l-1.2)     ■»    0.0 

IF      (1     .FU.     IS    ♦     I  )     60    10    41M) 

1     =     1-1 

r,o    r  n  h  1 7  o 

CONl  I  Vlr. 

IF      (IPI  IT     .LT.     31     Gil     10    h70 

WH  I  T  F      1  <'  •  1  1  0  0  ) 

FOwv.0)      (  111  I  lOXtlRMPTH     ARHAY      IN    Ft-0pl//| 

no    i 101    i    =    i ,   nc 

WHITF       I". 110?)      PTh(1i1)i     PTH(lt?)>     PlH(Ii3),     PlH(l,4) 
F"UHMA1     (2X,*mriaTf     IS     >F10.s.)h    A     ,  t  6  .  2  ,  I  OX  ,  mhS  T  APT     a     ,F10.?, 
IfthF'NlJ    =     >Elu.:>) 
»N !  1 t       111  I  04 ) 
FORMAT     ( IMOi lOXi IshKFTUNn    To    CpntS//) 

rum  i  nul 

Rt   1  IIHN 

FWi 


—  rnPT 

116 

FOPT 

1ST 

fopt 

158 

FopT 

\** 

FOPT 

lftO 

FOPT 

lftl 

FOPT 

1ft? 

FOPT 

Ift3 

FOPT 

1ft* 

FOPT 

lft5 

FOPT 

lftft 

FoPT 

lft7 

FOPT 

lftB 

FOPT 

lftO 

FOPT 

1T0 

FOPT 

171 

FOPT 

17? 

COSFOPT 

1T3 

FOPT 

17* 

FOPT 

175 

FOPT 

17ft 

FOPT 

177 

FOPT 

17* 

FOPT 

179 

FOPT 

l«0 

FOPT 

1*1 

FOPT 

1«? 

FOPT 

1«3 

FOPT 

)R4 

FOPT 

1*5 

FOPT 

lRft 

FnPT 

1*7 

FnPT 

1  HR 

FnPT 

1R9 

FOPT 

1^0 

FOPT 

191 

FOPT 

1«? 

FOPT 

103 

FnPT 

IP* 

FOPT 

|QS 

FOPT 

19ft 

OX,  FnPT 

197 

Fnp  t 

l^R 

FOPT 

199 

FOPT 

?no 

Fool 

?01 

FnPT 

?o? 

FOPT 

?nl 

SUHRhOuI'SM     I  E  NG7H. 
ll  U  '770 


IINUSF  D     rilMR  Il.fc   )     SH,\tf 
03 f 000 
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SLMhOIJTINE  ICALC  <A.R,n,C,H,AS,IS,HP,HA,l  EF.lTOT)                  ICAL 

C  CALCULATES  MOMENT  OF  INERTIA  FOr  A  PARTICULAR  SECTION  -   ICAL 

000016  REAL  ITOT,  LEF,  IPARTi  IS  ICAL 

C  CHECKS  FOH  COVER  PLATES  JCAL 

C  CALCULATES  I  FOR  NO  COVER  PLATES  —  ICAL 

000016  IF  IA  ,<iT,  ,000)  .OR,  0  ,GT,  ,00011  flO  TO  |A                        ICAL 

000027  IToi  •  IS  ICAL 

000030.  LEF  «  M^  ICAL 

000032  RE  r»)hN              »  ICAL 

C  »ETS  INITIAL  ZERO  VALUES  ICAL 

0Q0033  10  AL*  ■  0.0  -     ICAL 

00003*  A.M  ■  0.0  ICAL 

00003*  DL  »  0,0  ICAL 

000035  OU  •  0,0  ICAL 
C  CALCULATES  VALUES  FOR  THE  UPPER  PLATE  OR  SLAB                   ICAL 

000036  IF  (A  .LT.  .0001)  GO  TO  II  ICAL 
0000*1  AUt*  •  A  •  B  ICAL 
0000*2  DJ  *    MA  ♦  IH  ♦  R)/?.  ICAL 

C  CALCULATE  VALUES  FOR  THE  l.OWER  PLATE  ICAL 

0000*6  11  IF  lU  .LT.  .0001)  GO  TO  12  -     ICAL   2 

000051  ALP  ■  0  •  C  ICAL   2 

000052  DL  »  (H  ♦  C)/?.  ICAL   2 
00OOSS  12  ATOI  ■  AS  ♦  AUP  ♦  ALP  -   ICAL   2 

C  CALCULATE  DISTANCE  FROM  CFNTHOm  OF  STEEl  SECTION  TO  THE        ICAL   2 

c  ctNTKoio  of  the  overall  sfction,  also  caiculate  i  For  section,  ical  ? 

000060  DI S  I  •  (AUP  •  Dl)  -  ALP  •  DLI/ATOT  ICAL   2 

00006*  LEF  *  HP  ♦  UIST  ICAL   2 

000066  IPahT  .  |S  .  AUP.H»8/)2,  ♦  AlP*C»C/)?,  ICAL   2fl 

000076  ITOI  o  IPaRI  ♦  AUP'Ou"?  ♦  AlP»0L»«?  -  ATOT»MST»»2                ICAL   29 

000111  REfJhN  ICAL   30 

000112  ENO  ICAL   31 

SUHPROGfiAM  LENGTH 
0001*2 

UNUSED  COMPILER  SP«CE 
0*2500 


HEAL  FUNCTION  ILTNT  < PT , RO . COOR , I SW ,nP ,L FNflTH , SnBLEN) 
C         INTERPOLATES  REHfEN  POINTS  ON  ThF  INFIUfmCF  LlNE 

P0001?  OHtNSluN  RUIH1),  COoRIRlli  SUHLEN(RO) 

00001?  RtAu  LEHGIH 

000012  INlf.GFH  PI 

000012  IF  IIS*  .EU.  ?)  r.O  TO  5?0 

C         CALCUiATES  VMUE  OF  INFLUENCE  LlNF  TO  Thf  RIOhT  OF  POINT  PT 

00001*  SP  »  COOR(PT)  ♦  DP 

C         *tTUHNS  A  7fRO  IF  THE  POINT  IS  OFF  THF  OlODEd 


U  I  N    1 

II.  IN     I 

KIN   3 

|L  IN     « 

IL  IN    S 

1 L  I N    6 

IL  IN    7 

ILIN    A 

U  IN    9 

1  6  1 


00001*.  IF  (SP  .or.  LENGTH)  RO  TO  Sin 

0000??  GO  10  51? 

0000??  510  ILI'O  =  0.0 
0000?3  RFfoHN 

C         'JETEKMINES  Thf  SjaFLFME^T  LOCATION 

0000?!,  5'?  J  =  M  ♦  1 

0000?is  513  IF  icnow(.l)  .Gf  .  SP)  GO  TO  515 
00003?  j  =  J  ♦  1 

P00033  GO  1°  513 

00003".  515  K  =  J  -  I 

C         IMEKPOLATION  FORMULA 
000036  ILI-il  =  KUIK)  ♦  (HO(  j)  -oo(K)  )  «  ISP-COOPi' ) 1 /Sn U  f NIK) 

00005?  RM'iHN 

C         CALCULATES  v  ft  |  u  f"  OF  TNFLMFNCF  LINF  TO  Iwr  IFFT  OF  PojnT  PT 

00005?  5">0  S^  =  COOK(PT)  -  OP 

C         -<t  TUi'NS  A  /fmo  IF  THF  POINT  IS  OFF  Tht  cicOF) 
000054  IF  i',p  .L.T,  0.0)  RD  TO  5?? 

000056  GO  10  5?5 

000057  5??  ILM1  =  n,0 

000060  prf.iHN 

r  ilTFKMINES    Thf    S.IRFLFMFNT    LOCATION 

000061  5'5  J  =  PI  -     I 

000063  5?7  IF  ICOOH(J)  .U.  SP)  GO  TO  5?6 
000067  J  =  J  -  1 

000070  GO  u  5?7 

000071  5?6  K  =  J  ♦  1 
000073  J  =  * 
000071.  K.  =  J  -  I 

C         1 1  If  <POl_AT  ION  FORMULA 
000075  It  I  *1  =  R0(K)  ♦  (HO(J)-RO(K)  J  »  I  SP-fOOR  (K )  ) /Siml.EN  ( K ) 

O00I11  RFf.lfM 

00011?  END 


ILIN 

10 

ILIN 

1  1 

ILIN 

1? 

ILIN 

13 

ILIN 

1* 

ILIN 

15 

ILIN 

16 

ILIN 

17 

ILIN 

Ifl 

ILIN 

1" 

ILIN 

?0 

ILIN 

?1 

ILIN 

?? 

ILIN 

?3 

ILIN 

?4 

ILIN 

?5 

ILIN 

?ft 

ILIN 

?7 

ILIN 

?H 

ILIN 

?4 

an 

30 

1 1. 1  N 

31 

ILIN 

3? 

ILIN 

33 

|LtN 

34 

ILIN 

35 

ILIN 

3ft 

ILIN 

37 

ILIN 

3B 

ILIN 

39 

ILIN 

4(1 

ILIN 

4) 

SUBPHnGRAM    LtNGTH 
000133 

UNDSFO     CUMPII   ER     SP4C-F 
04?600 


S   HhIiUT  I  NF     ILF'KOP      (Rri.  S.  IRLF  N,NA  .KOP,  /VNFG.  "  P<">S  ,  M«XO|_  .MlNOl    ) 

OOOn  1  3  DIMENSION    KUIK1).     S-fRlFMlRO),     KOPI1M 

C  LFTMihlNFS    TmF     PROPFRTIF*    OF    A    r.TVFN    tmfiI'F'CE    |   TNF 

000013  MA  *   iL     =     1 

000013  mjxil  x  I 

00001 5  ANF  .  =  1,0 

000016  APi)->  =  n.O 

0  0  0  0  1  >,  0M4*  =  -.000] 

000017  0M  I  •«  =  .i001 
0000?!  DO  ->70  I  =  1,  10 
0000?3      5  '0  K1"l  I  )  =  n 

C  IfTEHt'I'NE  THF  CR3SSOVFP  POINTS 

0000?6  KOP ||)=1 

O0O0?7  J  =  r 


ILPH 

1 

ILPH 

? 

ILPu 

3 

ILPH 

A 

ILPH 

5 

ILPH 

6 

ILP-i 

7 

ILPH 

H 

ILPH 

4 

ILPH 

in 

ILPH 

1  1 

11  PH 

1? 

ILPH 

1  3 

ILPH 

14 
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oooo3ft 
00003? 
000034 
0000*3 
0000*=, 
0000<.7 
000047 
"00051 
00005ft 

ooooe3 

0000f<5 
000067 
00007? 
000074 

00007=, 
00007ft 
00010? 
000105 
000105 
000107 
0001 11 
0001 1? 
0001  14 
00011ft 

0001?1 

oool?? 

P00)?4 
P001?S 
O00l?7 
000131 
000133 
O0013S 
000137 
000145 
OOOlhO 
0  0  0  1  ft  1 
OOdftT 
OOOlhS 
O00173 
000?10 

oon?i  i 
ooo??o 

000??3 

ooo??s 

000?33 
000?33 
000?34 


5M 


550 


555 


554 


5*<? 


^3 


i;ft4 
5ft? 


5F0 
5*1 


NE 

D') 

IF 

KO 

J 

GO 

IF 

IF 

IF 

KO 

J 

CO 

KO 
NK 

0" 
IF 
IF 
60 
O'l 
MA 
GO 
OM 
MI 
CO 

Of) 
IF 
IA 
IH 

A 

no 

IF 

IF 

IF 

A 

GO 

IF 

IF 

IF 

A 

GO 

A 

CO 

IF 

IF 

Ci 

PF 

Ef 


a  NA 
-.SO  1 
(ARS< 
-MJ)  s 
J  ♦ 
10  55 
I  I  .F 
(H0(1 

I  ^H5  ( 

»»<  J)       D 
J      ♦ 

t I IMUF 

•M  J)  = 
M  =  J 
•  lETE 
-.5?  I 
IKI(  I 
[Ml  (I 
fO    5^ 

A*     =     . 

K.Jl      = 

10  5'. 

I    i       =       H 
I'll        I 
-J  I  I 'Mill-. 

IF  TF" 
->ftO  .K 

IK  .F 
=  IU)R 
a    KOR 

*  o.o 
-.ft?   I 

II  ."1 

I  I  .F 
iahSI 

=     bllHL 

11  S< 

I  1     .1. 

HO 

IftHS  ( 
x    A«SU 

10    5  ft 
i    A     ♦ 
II 1NUF 

I  A     .(, 

(A     .L 
Nl  INIIF 
I    it-fj 
>) 


-     1 

*     <?.     NF 
RO(I))     .r,T.     10.0»»l-ft>)     GO    TO    551 

1 
1 
o 

0.     1  >    (.0    TO    550 
)»RO(I-l)     .GT.    n.o)    (50    TO    ^5n 
HO(I-])>     .LT.     lo.0»»<-6))     GO    TO    55n 

I 
1 

MA 

MINE    MAMMIIM    AND    MINIMUM    OPOINATF     I   "NATIONS 

=     1  .     Fl  A 
)     .GT.    OMAX)     GO    TO    555 
I     .LT.    OmIM)     GO    TO    554 

(MI) 

I 

? 

(1(1) 

I 

MINE  ThF  POSITIVE  ANO  NEGATIVE  ARFA5  IJNnfH  THE  LINE 
=  1 ,  NKOP 

0.  NKf'P)  (.0  TO  R6I 
<K) 

("♦II  -  1 

-      IA.  IP 
F .  I A I  GO  TO  5ft3 
').  ll  GO  TO  c,63 

"("■(I))  .1  T.  )0.0*»(-ft))  CO  TO  c,fti 

FMI-])<>POII)«AHS(c!0(I))/<?.»(AP5(Rn,T))4ARS(RO(l-l)>)) 
4 

F  .   IH)  GO  TO  c;f,4 
,EQ.  MA)  GO  TO  c.ftft 

K0(I*1))  ,lT.  lo,0»»<-6))  GO  TO  Cft4 

RLENI I ) »RO ( I )«ABS(Rn< I) >/ <?.»(»R5 (RO( 1) > ♦ARSfRol T«l ) > ) ) 
? 
(r(0(ii  ♦  »0  ( i  ♦  l ) )  •  SUSLFN)  T )  /?. " 

1.  0.0)     aRo5    =    A»05    .    A 
t.     0.0)      ANFG     b     ANFG     ♦     A 


ILPR 

15 

ILPR 

16 

ILPR 

17 

ILPR 

lfl 

ILPR 

19 

ILPR 

?0 

ILPR 

?1 

ILPR 

??. 

ILPR 

?3 

ILPR 

2* 

ILPR 

?5 

ILPR 

?6 

ILPR 

?r 

ILPR 

?R 

ILPR 

?9 

ILPR 

30 

ILPR 

31 

ILPR 

3? 

ILPR 

33 

ILPR 

3* 

ILPR 

35 

UPR 

36 

ILPR 

37 

ILPR 

3B 

ILPR 

39 

ILPR 

40 

ILPR 

41 

ILPR 

42 

ILPR 

43 

ILPR 

44 

ILPR 

45 

ILPR 

46 

ILPR 

47 

ILPR 

4H 

ILPR 

49 

ILPR 

50 

ILPR 

51 

ILPR 

5?. 

ILPR 

53 

1LP.J 

54 

JLPw 

55 

ILPR 

5ft 

ILPR 

=.7 

ILPW 

5R 

II  PR 

59 

ILPH 

ftO 

ILPR 

ftl 

ILPW 

ft? 

II  PR 

ft3 

SUHPRO  'RAM     LF  NbTrl 
00O?77 


UNIJSFO     TO-'PII   FH     S^nl  F 
04??U0 
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RC»i.    FUNCTION    Impact     (Type  .("P.  IOSP,*  Si«ci.r  ran. COOS) 
C  oITEWminES    Thf     Ixp»CT    FACTUM 

00001?  D!  ""'NSInrj    IOSP(^).    SpANii,).    rOCmifln 

00001?  INfMiFH    TYPE 

00001?  MN  a  NS  -     1 

000013  U  ""'N  •E,J»  11  c.o  rn  7nno 

OOOOIS  IF  (TrPf  .NE.  (.hpFAC)  (40  TO  600 

C  IMPACT  DISTANCE  Fort  ThF  RFACTIOn  TNFI'iFMrr  LINES 

000017  IF  INP  «Nt.  1)  Rn  Tf)  <i01 

0000?0  OlN'>  =  SPANd  I 

000021  GO  (0  650 

0000??  f.0 1  IF  INP  .NE.  NM  00  To  6m? 
0000?4  D  I  M  •>  »  SPAN(NS-l) 

0000?6  GO  lO  6SI' 

0000?*  602  D I  ^>i  -    SPAN(NP)  *  SPaN(NP-I) 
00003?  GO  10  650 

00003?  600  If  (TYPF  .NE.  4HMOMT)  Go  TO  610 

C  It-PACT  DISTANCE  FOR  THE  MOMENT  IMFLUFmCF  i  TNfS 

000034  DO  ill  I  e  2.  NS 

000035  IF  IIOSP(I)  .LT.  NO)  (JO  TO  611 

000040  IF  ll'fAf)  .LT,  O.il  Go  To  61? 

000041  DI*»  =  SPAN! I-i ) 

000043  GO  10  65(J 

000044  61?  IF  !&P  .LT.  ( IOSP(  I-l ) *TOSP( I  I ) /?)  oo  To  *?c 
000051  DIi.i  =  <5PAN(1-J1  •  SpAN(I)»/2.0 

000056  GO  10  650 

000056  6?5  01  -J  .  =  (5PANII-1)  ♦  SPAN ( I -? ) I /2 . 0 

000063  GO  TO  6Sp 

000063  611  COmiIMUE 

000066  610  IF  llVPF  . NE .  4HSHFR)  GO  TO  6?0 

C,         IMPACT  DISTANCE  FOR  THE  SHEAR  IMFLUFNCF  i  tNFS 

000070  DO  »?1  1  =  2.  NS 

000071  IF  (IOSP(I)  .IT.  NO  GO  TO  6?1 

000074  IT  =  lOSP(I) 

000075  J.I  =  IQSPlI-ll 

000077  DI M'i  =  COOR(II)  -  COOR(NP) 

000104  DO*(.>  =  COOfl(NP)  -  COOR(JJ) 

00C110  IF  iDONo  .GT.  DING)  OlNf,  ■  DONG 

000113  GO  10  650 

000114  6?1  COMIlNUF 

C  IMPACT  DISTANCE  FOR  THE  DEFLECTION  TNFUIF-CF  LlNFS 

000117  6?0  00  h3o  I  =  2.  NS 

0001?1  IF  <NP  .GT.  lOSP(I))  GO  TO  630 

0001?5  D I  Mi  =  SPAN(l-) ) 

000l?6  60  10  6S() 

000127  630  COMflNllF 
00013?  GO  10  650 

000132     7000  DI  M  i  =  SPAN! 1 ) 

C         lf-PACT  FACTOR  CALCULATION 

000134  6C0  HT  »  50. /(  (DING/1?. )  ♦  1?S.) 
000137  IF  IhT  .GT.  .3)  hT  r  .3 

000143  HT  a  HT  ♦  1.0 

000145  IMVaCT  =  HT 

000146  REIJhN 
000146  END 


I"PT 

1 

I'<PT 

? 

IMpT 

3 

IMP! 

4 

I"PT 

5 

MPT 

6 

IMPT 

7 

IMPT 

H 

IMPT 

9 

MPT 

10 

IMPT 

11 

IMPT 

1? 

IMPT 

13 

IMPT 

14 

IMPT 

15 

IMPT 

16 

IMPT 

17 

IMPT 

18 

IMPT 

19 

IMPT 

20 

IUPT 

21 

IMPT 

22 

I"PT 

23 

I  <PT 

2* 

IMPT 

25 

riPT 

26 

IMPT 

27 

IMPT 

2H 

IMPT 

29 

MPT 

30 

IMPT 

31 

IMPT 

32 

IMPT 

33 

IMPT 

3* 

IMPT 

35 

IMPT 

36 

I  "PT 

37 

IMPT 

3H 

IMPT 

39 

IMPT 

40 

IMPT 

41 

I'iPT 

42 

IMPT 

43 

IMPT 

44 

IMPT 

45 

IMPT 

46 

IMPT 

47 

IMPT 

4fl 

IMPT 

49 

IMPT 

50 

IMPT 

51 

IMPT 

52 

I  1PT 

53 

IMPT 

54 

IMPT 

55 

1  6 


SUHPHlriRAM    |_FM.'TH 

ooo?n 

UNUSED    rU^PIlFK    Sr'nCF 
042*00 


SUljrfOU 

C  ->'U 

oooon  coiicn 

000011  CCMniON 

OOOOll  DUi-.NS 

C  CHF 

oooon  if    up 

000017  IF     lIC 

000021  RE'  "HN 

000021             7003  If     (*T 

0000?*  RE  MhN 

C  »T0 

00  0  025            70 r  4  DO     /GO 

0000?7  DO    "SI 

000030  PLT-iAV 

000036            6919  PLWiAV 

0000*3            70H0  IS4vt ( 

0000*7  CHI i    = 

000051  WMI <    = 

000051  00     '00 

000053  00     100 

00005*  00     (00 

000055            70''1  SAVFII 

00007IS  RE  r  iRN 

000077  EN) 


IT1NF     KFEPfNE .NOS.PLATh.PLWID.COVPI   ,CMTm,UmIN,  IPTft) 

iRFS    ThF    ArrFOTAHLF    GIRDER 
)/f  OUR/wTOT.TC><;T 


l/l-  OUR/^TOT.TCXiT 

l/A/lS«VE(«0>.SAVF(l2»6.2).PLT5«V<nn,?i.P|.wSAv(«0.?> 
ION    NU5(flp)  ,PLATH(R0i?l  •PLWin(OTO>?)  .rr"'P|   (12, S, 2) 
C*     IF    THF     r,I«TFR    15    *ETTFR    Than    Thf    PorVlnuS    RTRnFo 
Th    ,QT.     .*•)     RO    TO    7003 


OST     ,LT.    CmIN)     till    JO    700* 
OT    AT,    hMTNl     r,0    TO    700* 

RE     THE    ACCFPTAflLF    RIROFR 

0  I     =     1  .     NF 
9     J     =     1  .     ? 
(I.J)      e     PL  ATHl  J  ,  J) 

1  1 . J  )     =    f-  L  W  I  0  (  I  ,  J  ) 
I)      *     NOMI) 

TCOST 
*TOT 
1    I    »    li    1? 

1  J  "  1  .  ft 
1  K  a  1,  ? 
, J.KI    a    COVPL ( I i J.K) 


SUBPROGRAM  LFN^TH 

000120 


UNUSED  (-OHPII  ER  S^aCE 
0*2500 


K.FFP 

1 

KFFP 

2 

KEEP 

3 

KEEP 

* 

KFFP 

5 

XFFP 

6 

KEEP 

7 

KFFP 

R 

KEEP 

9 

KFFP 

10 

KEEP 

11 

KEFP 

12 

KEEP 

13 

KEEP 

1* 

KEEP 

15 

KEEP 

16 

KEEP 

17 

KEEP 

18 

KEEP 

19 

KEEP 

20 

KEEP 

21 

KFEP 

22 

KEEP 

23 

KEEP 

2* 

KEEP 

25 

000021 

000021 
000021 
000021 
000021 

000023 


.LFMG 


TH, COOT, NP. TYPE. LOAD 
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